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INTRODUCTION 


The  Rotorcraft  Flight  Simulation  Program  C31  is  a modular 
computer  program  capable  of  performing  performance,  stability, 
and  rotor  loads  solutions.  This  complex  and  comprehensive 
program  may  be  operated  in  a variety  of  ways.  Data  input  and 
selection  of  user  options  determine  the  degree  of  interaction 
between  the  different  analyses  offered  by  the  program,  and 
determine  the  complexity  of  each  of  the  major  sections  of  the 
program.  The  C81  essentially  utilizes  the  modular  concept  of 
programming  the  different  major  analyses  comprising  a global 
program  package.  The  program  is  described  in  some  derail  in 
References  1 through  3. 

The  C81  program  was  developed  by  the  Bell  Helicopter  Company, 
and  has  been  made  available  to  the  rotorcraft  industry  by  the 
U.  3.  Army.  Although  used  and  validated  mostly  with  helicop- 
ters with  teetering  rotor  systems,  the  program  is  also  intended 
for  modeling  helicopters  with  articulated  or  hingeless  rotors. 
At  the  request  of  the  U.  S.  Army  Air  Mobility  xResearch  a.nd 
Development  Laboratory,  Lustis  Directorate,  this  study  was 
conducted  by  Sikorsky  Aircraft  to  validate  the  C81  program 
for  use  with  articulated  rotor  helicopters. 

This  report  describes  the  comparison  of  the  C81  program  with 
test  data.  Correlation  results  are  presented  and  discussed 
for  the  H-53  (S-65)  and  S-67  aircraft.  The  usefulr.ess  of  the 
C81  program  as  an  articulated  rotor  helicopter  design  tool  is 
evaluated  based  on  the  correlation  data  gathered  in  the  study. 
This  report  also  provides  constructive  criticism  of  those 
areas  of  the  program  needing  improvement  to  increase  accuracy. 

Both  the  Sikorsky  K-53  and  the  Sikorsky  S-67  aircraft  represent 
single-rotor  helicopters  with  anti-torque  tail  rotors.  The 
main  rotors  of  these  aircraft  are  fully  articulated,  employing 
blade  flapping  and  lagging  (hunting)  hinges.  The  twin  turbine 
engine  H-53  is  basically  a transport  helicopter,  although  it 
has  also  been  used  for  a number  of  other  purposes.  Many 
production  versions  of  the  H-53  exist,  but  the  basic  charac- 
teristics of  these  versions  are  quite  similar.  Figure  1 
is  a drawing  of  the  CH-5  3A.  The  tv/in  turbine  engine  S-67 
demonstrator  aircraft  was  designed  as  a high-speed  derivative 
of  the  Sikorsky  S-61  (SH-3D)  helicopter.  Wings  could  be 
attached  to  the  sponsons  for  additional  lift  and  attachment 
points  for  armament.  Figure  2 is  a drawing  of  tlie  S-67. 
Principal  dimensions  and  general  data  for  the  H-53  and  S-67 
aircraft  are  presented  in  Table  I. 
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PROGRAM  IMPLEMENTATION 


The  Eustis  Directorate  provided  Sikorsky  Aircraft  with  a copy 
of  the  300K  (300,000)  bytes  of  core  storage  rendition  of  the 
Rotorcraft  Simulation  Program  C81,  Version  AGAJ74.  This 
version  of  the  program  was  adapted  for  use  on  the  contractor's 
IBM  370/155  computer  system.  The  300K  core  version  of  the  C81 
program  differs  from  the  full  size  600K  core  version  of  the 
program  primarily  in  the  number  of  input  data  tables  that  may 
be  loaded  into  the  program.  Rotor-induced  velocity  of  the 
rotor  wake  acting  at  aerodynamic  surfaces  cannot  be  loaded  into 
the  smaller  300K  version  of  the  program.  These  phenomena  are 
described  instead  by  analytic  expressions  in  the  version  of  the 
C81  used  in  this  correlation  study.  Whereas  the  BOOK  version 
of  the  C81  provides  for  up  to  five  airfoil  data  tables  to  be 
loaded  into  the  program  to  describe  the  aerodynamic  character- 
istics of  rotor  blades  or  aerodynamic  surfaces,  the  300K  version 
of  the  C81  has  room  for  only  two  such  tables. 

No  major  problems  were  encountered  in  converting  the  IBM  360-65 
copy  of  the  program  for  operation  on  the  IBM  370/155  computer 
system.  The  program  was  initially  compiled  using  a G2  compiler. 
Later,  in  an  effort  to  reduce  the  running  time  of  the  program, 
it  was  re-compiled  using  the  H2  (H  - extended)  compiler.  (The 
running  time  requirement  is  discussed  in  the  General  Program 
Characteristics  section  of  this  report.) 

A number  of  check  cases  were  provided  by  USAAMRDL  to  determine 
if  the  C81  program  had  been  converted  properly  to  the  IBM 
370/155  system.  The  check  cases  included:  (1)  a quasi-static 

trim  followed  by  a quasi-static  maneuver  of  approximately  two 
seconds  duration  for  the  AH-IG;  (2)  a quasi-static  trim  followed 
by  a time-variant  trim  and  maneuver  of  approximately  four 
seconds  duration  for  the  AH-IG  with  an  automatic-pilot  ''eedback; 
and  (3)  a quasi-static  trim  and  stability  solution  of  a modified 
UH-lD/H  deck  which  analytically  simulated  an  articulated  rotor 
helicopter.  The  input  data  for  these  cases  were  loaded  into 
the  IBM  370  version  of  the  C81  and  the  program  was  executed. 

The  solutions  were  compared  with  the  check  case  solutions 
provided  by  the  Eustis  Directorate.  The  correlation  of  these 
three  check  cases  was  discussed  in  detail  with  the  contract 
technical  representative  at  E’ort  Eustis;  the  printout  for 
check  case  1 was  delivered  to  the  Eustis  Directorate  as  required 
by  contract.  With  the  exception  of  the  portion  of  the  time 
history  mcineuver  response  which  followed  the  activation  of  the 
Automatic  Pilot  routine,  all  portions  and  options  of  the  C81 
program  investigated  by  the  check  cases  at  the  contractor's 
facilities  yielded  the  same  results  as  those  provided  by  the 
Eustis  Directorate.  The  only  inconsistency  involved  the  Auto- 
matic Pilot  routine  and  was  explained  as  an  error  or  modifica- 
tion in  the  routine.  But, since  the  Automatic  Pilot  routine 


would  not  be  needed  for  the  planned  correlation  study,  it  was 
decided  that  the  C81  program  had  been  made  properly  operational 
at  the  contractor’s  facilities. 

About  midway  through  the  C81  prediction/correlation  effort,  a 
CALCOMP  plotting  system  was  installed  at  the  contractor's 
facilities  and  was  interfaced  with  the  C81  program  which 
possesses  a CALCOMP  plotting  capability.  The  CALCOMP  routine 
was  used  to  plot  some  of  the  time  histories  that  were  run  in 
the  correlation  effort.  An  evaluation  of  the  C81  CALCOMP 
routine  can  be  found  in  the  Discussion  of  General  Program 
Characteristics  section  of  this  report. 


PROGRAM  INPUT  DATA  PREPARATION 


Rotor  blade  airfoil  data,  aerodynamic  data,  blade  bending  mode 
shape  data,  and  all  other  necessary  aircraft  characteristic 
data  required  as  input  to  the  C81  program  were  prepared  for 
the  H-53  (S-65)  and  S-67  single-rotor  helicopters.  The  H-53 
and  S-67  input  data  for  the  C81  program  were  drawn  from  Refer- 
ences 4 through  9.  Some  of  the  basic  data  from  these  references 
were  manipulated  as  required  to  present  them  in  a form  compa- 
tible with  the  C91  program.  The  digital  computer  program 
AS812A,  written  in  the  PL/1  computer  language,  was  obtained 
from  the  Eustis  Directorate  and  used  to  convert  fuselage  wind 
tunnel  data  into  aerodynamic  input  data  for  the  C81  for  the 
H-53  and  S-67  aircraft.  An  existing  in-house  Sikorsky  digital 
computer  program  was  used  to  calculate  the  blade  normal  mode 
shapes.  These  data,  in  turn,  were  manipulated  by  hand  to 
prepare  them  for  the  C81  input  format.  A short  FORTRAN  pro- 
gram was  written  to  convert  existing  airfoil  data  from  th^ 
Sikorsky  input  format  to  an  input  format  compatible  with  the 
C81.  Appropriate  airfoil  data  were  then  converted  for  the 
correlation  effort.  The  Sikorsky  airfoil  data  conve..ted  for 
use  with  the  C81  program  consisted  of  data  which  were  smoothed 
to  best  fit  a nximber  of  sets  of  raw  data  obtained  from  many 
wind  tunnel  tests  of  each  airfoil  section  under  consideration. 
Comments  regarding  the  C81  input  format  and  problems  en- 
countered while  preparing  data  for  input  into  the  C81  program 
are  presented  in  the  Discussion  of  General  Program  Character- 
istics section  and  other  related  sections. 

Tables  II  through  VII  are  listings  from  the  C31  program  of  the 
basic  inputs,  airfoil  data,  and  blade  aeroelastic  data  for  the 
H-53  and  S-67  aircraft  which  were  used  in  the  majority  of  the 
C81  runs  conducted  as  part  of  the  correlation  effort.  When  any 
input  datum  is  related  to  a problem  or  an  inadequacy  in  the  C81 
program,  the  datum  is  described  in  sufficient  detail  in  the 
appropriate  sections  of  the  report. 
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■ ll-^TllOJ  OF  COHRKI^VnOO  A.ND  DFSC RIPTIO.-I  OF  TFST  DATA 

The  validation  of  the  C31  program  was  accomplished  through 
comparison  of  test  data  with  the  data  calculated  by  the  COl 
rotorcraft  analysis.  For  each  correlation  case  investigated, 
helicopter  speed,  weight  loading,  and  control  inputs  were 
loaded  into  the  C81  to  best  duplicate  the  conditions  flown  in 
the  test  program.  The  C81  and  test  data  correlation  results 
were  drawn  on  to  evaluate  the  accuracy  of  the  C81  program  in 
modeling  articulated  helicopter  static  and  dynamic  flight 
cliaracteristics . 

Flight  test  and  whirl  stand  test  data  used  in  the  correlation 
effort  were  selected  based  on  confidence  in  the  accuracy  of 
the  data  as  well  as  the  ease  of  accessibility.  The  test  data 
represented  cases  chosen  to  exercise  the  three  major  ['ortions 
of  the  C81:  trim,  stability  and  control,  and  maneuver.  Where 

available,  results  from  other  digital  rotorcraft  analyses  were 
also  used  for  comparison  with  the  C81  program. 

The  major  prediction  efforts  conducted  under  this  study  are 
more  specifically  defined  by  five  types  of  data  comparisons: 

(1)  Trim 

(2)  Performance 

(3)  Stability  (roots  and  derivatives' 

(4)  'I’ime  history  (helicopter  response  characteristics) 

(5)  Rotor  loads 

Table  VIII  lists  all  of  the  conditions  under  these  five  areas 
that  were  successfully  executed  by  the  C81  program,  producing 
data  for  the  correlation  study.  The  reports  from  which  the 
data  were  selected,  against  which  the  C81  is  compared,  are 
listed  in  the  reference  section.  Pertinent  test  and  analytic 
data  from  these  sources  and  the  corresponding  C81  data  appear 
in  the  following  sections  of  this  report. 

The  test  data  used  in  the  correlation  effort  were  taken  from 
Reference  7 and  from  References  9 through  16.  The  trim  data 
from  References  9 and  10,  and  the  unpublished  d-67  trim  data 
are  raw  flight  test  data.  The  unpublished  S-67  performance 
data  and  the  forward  speed  performance  data  from  Reference  12 
also  are  raw  flight  test  data.  The  hover  performance  data 
from  Reference  12  are  test  data  which  have  been  analytically 
fitted  and  extrapolated  from  raw  flight  test  data.  The  main 
rotor  and  tail  rotor  whirl  stand  data  are  refined  from  raw  test 
stand  data  to  reflect  standard  sea  level  conditions.  The  climb 
and  autorotation  test  data  from  Reference  14  are  refined  from 
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raw  flight  test  data  to  reflect  consistency  in  atmospheric 
conditions.  The  root  locus  data  from  References  15  and  16 
hc.ve  been  processed  from  smoothed  time  history  flight  test 
data.  All  time  history  and  rotor  loads  data  from  References 
7,  11,  14,  and  16  have  been  filtered  to  remove  high  frequency 
noise  content  from  the  flight  test  data.  All  test  data  selected 
here  for  correlation  purposes  are  of  high  accuracy  within  each 
specific  data  type  (trim,  performance,  response,  and  rotor 
loads) . 

Contract  schedule  and  budget  restrictions  dictated  the  number 
of  C81  runs  made  during  the  study.  (Approximately  260  explora- 
tory and  data  runs  were  made.)  All  C81  program  updates  and 
modifications  specified  by  the  program  developer,  through 
USAAMRDL,  for  the  AG.AJ74  version  were  incorporated  in  any  H-53 
and  S-67  data  runs.  Flight  conditions  that  were  attempted  but 
not  run  successfully  (and  therefore  not  listed  in  Table  VIII) 
are  discussed  in  the  following  sections.  Those  flight  condi- 
tions in  Table  VIII  that  were  run  with  only  moderate  success  or 
v’th  great  difficulty  are  also  elaborated  on  in  the  following 
sections . 

1 

, The  following  sections  of  the  report  group  the  results  and 

I evaluation  of  the  C81  comparative  study  according  to  the  major 

I analyses  conducted  by  the  C81  program.  Any  noteworthy  charac- 

teristics of  the  C81  experienced  while  running  or  analyzing 
particular  aspects  of  the  program  are  described  in  the  appro- 
priate sections.  Good  and  poor  characteristics  of  the  C81  that 
have  been  discovered  in  this  correlation  effort  are  discussed. 
The  activation  or  deactivation  of  levels  of  program  sophisti- 
, cation  or  the  selection  of  program  options  is  defined  or  ex- 

I plained  when  necessary.  Also,  a section  of  this  report  is 

t devoted  to  discussion  of  general  program  input,  operational, 

] and  output  characteristics  that  are  not  specifically  related 

j to  one  of  the  five  basic  types  of  program  analyses  under 

investigation. 


TRIM  CORRELATION  - RESULTS  AND  DISCUSSION 


The  general  trim  solution  of  the  C81  program  is  of  great 
interest,  because  practically  every  type  of  analysis  offered 
by  the  program  begins  by  calculating  the  trim  conditions  oi 
the  rotorcraft.  In  computing  data  cases  for  the  articulated 
rotor  correlation,  a few  problems  were  experienced  perteining 
to  trim  iteration  solutions.  These  occasional]''  experienced 
difficulties  are  discussed  f,vSt. 

Although  the  version  of  the  C81  program  used  in  t is  validation 
study  contained  its  own  table  of  0012  airfoil  data,  another 
table  was  prepared  in  the  C81  input  format  based  on  0012  air- 
foil data  existing  at  the  contractor's  faciJity.  Also  converted 
for  input  to  the  C81  program  was  a table  of  modified  OOil  air- 
foil data,  intended  for  use  in  modeling  the  H-53  aircraft. 

Often  initial  efforts  to  trim  the  H-S^,  modeled  by  the  C81, 
would  not  succeed,  even  following  as  many  as  "iO  trim  iterations. 
Study  of  the  trimming  routine  partial  derivative  ixatrix  revealed 
that  the  main  rotor  moments  were  oscillating  by  repeatedly 
identical  values,  above  the  trim  tolerance  loadeu  into  the 
program.  When  the  contractor's  version  of  0012  ctirfoil  data  was 
used,  occasionally  the  same  problem  would  occur.  When  the  C81 
version  of  0012  airfoil  data  was  substituted  in  an  otherwise 
identical  case,  trims  were  always  successfully  obtained,  even 
though  differences  in  the  two  sets  of  0012  data  were  apparently 
negligible.  In,  order  to  run  the  progrcm  as  efficiently  as 
possible,  and  because  H-53  flight  characteristics  were  shown 
not  to  be  highly  influenced  by  the  dilferenoss  in  0011  and  0012 
airfoil  data,  the  C81  version  of  the  airfoil  data  was  used 
exclusively  during  the  remainder  of  the  correlation  study. 
Although  it  was  not  determined  why  the  of  only  the  C81-sup- 
plied  airfoil  data  guaranteed  successful  corvergence  of  trim 
rotor  moments,  a possible  explanation  ma',  be  that  the  airfoil 
angle-of-attack  and  Mach  number  increments  selected  for  loading 
the  airfoil  data  table  at  high  angles  cf  attack  are  more  criti- 
cal than  originally  thought.  Tables  IV  and  V may  be  used  to 
compare  the  0012  airfoil  da  a supplied  within  the  C81  program 
to  the  0012  airfoil  data  prepared  for  the  C81  from  Sikorsky 
data  tables  for  this  airfoil.  The  C81-supplied  0012  data  tables 
were  used  in  modeling  the  S-67  by  the  C81  program  in  order  to 
avoid  re-experiencing  trimniing  difficulties  similar  to  those 
encountered  witli  the  H-5  3 simulation. 

In  rare  instances  it  was  not  possible  to  trim  the  helicopter 
while  fixing  either  the  roll  or  yaw  angle  at  a specific  value. 
This  situation  was  experienced  while  trying  to  trim  the  S-67  at 
a 1.5-deg  yaw  angle.  The  C81  would  not  trim  the  helicopter  for 
either  +1.5  deg  yaw,  yet  it  trimmed  successfully  at  a zero  yaw 
angle.  Side  force  and  yawing  moment  were  not  different  enough 
between  these  different  yaw  conditions  to  account  for  tlie 
experienced  trim  difficulties. 


17 


For  certain  flight  conditions,  the  C81  trim  routine  appeared 
to  be  very  sensitive  to  the  input  starting  values  of  control 
position,  particularly  the  pedal  setting.  For  some  runs,  the 
trim  routine  successively  modified  the  control  input  in  the 
^j"ong  direction,  without  being  able  to  correct  itself  within 
the  given  number  of  trim  iterations,  and  thus  failed  to  trim. 
This  situation  occurred  quite  frequently  when  the  C81  attempted 
coordinated  turns,  and  occasionally  occurred  when  the  C81 
attempted  various  other  flight  conditions. 

Sometimes  the  main  rotor  flapping  angles  a^  and  b^  were 

printed  out  as  asterisks  (*)  in  the  time-variant  trim  solution. 
Experience  running  the  program  indicated  that  this  occurred 
whenever  the  value  of  XIT(2)  {the  rotor  azimuth  increment  during 
the  time-variant  trim)  , was  not  small  enough  to  allov/  at  least 
ten  calculations  per  highest  blade  bending  mode  frequency.  This 
characteristic  was  observed  a number  of  times  even  though  the 
user's  manual  claims  that  the  C81  will  calculate  and  substitute 
the  proper  value  for  the  parameter  XIT(2)  if  an  improper  value 
is  loaded  into  the  program.  The  lack  of  a solution  for  a^  or 

bj^  in  a time-variant  trim  prevents  the  program  from  proceeding 

to  following  stacked  cases  (if  requested). 

In  isolated  cases,  when  at  least  five  bending  modes  were  input 
for  a quasi-static/time-variant  trim  search,  the  time-variant 
rotor  flapping  angles  printed  out  were  asterisks.  There  pro- 
bably was  not  enough  memory  storage  for  the  program  (particu- 
larly if  there  had  been  a high  numb-„r  of  rotor  blades)  , so  that 
some  values  were  overwritten,  tlius  -‘"■ntaminating  portions  of 
the  solution. 

Another  general  trim  problem  concv.t",ed  stacking  cases  using 
different  numbers  of  rotor  elastic  modes.  The  C81  would  not 
converge  on  a quasi-static  trim  if  the  preceding  case  was  a 
quasi-static  trim  for  a condition  using  a higher  number  of 
elastic  modes.  Most  likely,  some  residual  variables  in  the 
trim  solution  of  the  preceding  case  associated  with  the  blade 
mode  shapes  were  not  zeroed  out  before  the  C81  began  the 
unsuccessful  case. 

Turning  now  to  the  specific  trim  conditions  tliat  were  studied, 
the  buliv  of  the  C81  trim  data  cases  that  were  run  requested 
the  quasi-static  trim  option  and  were  run  with  only  two  main 
rotor  elastic  modes  (hinged  rigid  flanping  and  lagg  iig)  . Ini- 
tial C31  runs  that  were  made  to  investigate  the  effects  of 
five,  four,  or  two  rotor  normal  modes  showed  that  tJie  number 
of  modes  had  little  influence  on  trim  values  of  aircraft 
attitudes,  control  positions,  or  flapping.  Most  pertinent  trim 
parameters  are  held  constant  within  the  C8I  when  going  from  a 


quasi-static  to  a time-variant  trim  solution.  Investigatory 
runs  showed  negligible  variation  in  trim  parameters  between 
the  quasi-s tatic/tii:ie-variant  pairs  of  solutions.  Yawed  flow 
and  unsteady  aerodynamics  also  proved  to  have  negligible  effect. 

Figures  3,  4,  and  5 show  the  trim  pitch  attitude  control  angle 
versus  speed  for  the  CH-53D  modeled  by  the  C81,  compared  to 
flight  test  data  at  three  center-of -gravity  positions.  The 
collective  pitch  calculated  by  the  C81  is  typically  higher 
than  that  indicated  by  flight  test  data  throughout  the  speed 
range,  except  at  high  speeds,  where  it  is  lower.  Tne  difference 
in  tail  rotor  pitch  between  the  test  data  and  the  values  pre- 
dicted by  the  C81  reflects  the  inability  of  the  C81  to  model 
the  delta- 3 (flap/pitch)  coupling  of  the  tail  rotor  blades  of 
the  CH-53D.  The  delta-3  was  loaded  into  the  program  input,  but 
the  program  solution  was  identical  to  a case  run  using  no 
delta-3.  The  delta-3  causes  the  effective  blade  pitch  to  be 
lower  than  the  pitch  impressed  by  the  control  system.  The  test 
data  are  impressed  control  system  inputs?  the  values  calculated 
by  the  C81  are  effective  mitch  angles. 

Figure  6 shows  S-67  trim  data  versus  speed  for  mid  center- 
of  -gravity  loading,  while  Figure  7 shows  other  trim  parameters 
for  the  S-67  at  the  forward  center  of  gravity.  Both  figures 
contain  flight  test  and  C81  predicted  data.  The  fuselage 
angle-of-attack  correlation  in  Figure  6 illustrates  a C81 
program  deficiency  in  that  the  influence  of  the  rotor-induced 
velocity  on  the  fuselage  is  not  calculated.  Also,  the  lack  of 
such  a downwash  calculation  certainly  has  some  influence  on 
the  accuracy  of  prediction  of  the  required  thrust  and  main  rotor 
collective  blade  pitch. 

In  general,  the  data  of  Figures  6 and  7 indicate  that  the  predic- 
tion of  aircraft  attitude  and  cyclic  control  angles  by  the  C81 
becomes  less  accurate  at  increasing  forward  speeds.  The  program 
solution  of  roll  attitvide  in  particular  seems  to  be  overly 
sensitive.  Similar  to  the  CH-5  3D  le^^el-f light  correlation  of 
Figures  3,  4,  and  5,  the  S-67  level-flight  correlation  indicates 
that  the  C81  predicts  lower  collective  pitch  requirements  at 
high  speed  than  are  indicated  by  the  flight  test  data.  The 
rotor  aerodynamic  characteristics  at  high  local  blade  angles  of 
attack  may  not  be  calculated  with  sufficient  accuracy  by  the 
C81.  Unlike  the  CH-53D  correlation,  the  S-67  data  comparison 
shows  the  required  C81  trim  collective  pitch  at  low  and  moderate 
speeds  to  be  lower  than  the  te.st  values.  The  wing  of  the  S-67 
may  be  carrying  too  much  lift  in  these  C81  simulations.  Similar 
to  the  CH-53D  correlation,  the  S-67  tail  rotor  pitch  requirement 
predicted  by  the  C81  differs  somewhat  from  the  flight  test  data 
due  to  the  delta-3  hinge.  For  the  S-67,  another  contribution 
to  the  difference  in  tail  rotor  pitch  may  be  due  to  a less  than 
accurate  calculation  of  the  force  produced  by  the  highly  effec- 
tive vertical  tail  of  this  helicopter. 
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rhe  prediction  of  the  S-67  pitch  attitucie  and  lomjitudinal 
cyclic  control  between  60  and  100  kt  forward  speed,  seen  in 
Figure  7 describes  an  unusual  trend  in  the  values  predicted  by 
the  C81  program.  The  jump  in  these  data  might  suggest  a severe 
discontinuity  in  either  the  fuselage  aerodynamic  data  or  the 
airfoil  data  inputs;  however,  no  such  discontinuity  was  dis- 
covered upon  examination  of  these  data  tables,  A possible 
explanation  for  this  isolated  trim  characteristic  may  lie  in  the 
aerodynamic  calculations  of  the  wing  or  horizontal  tail  within 
the  C81  program. 

Inconsistencies  between  C81  and  test  trim  parameters  are 
influenced  by  the  aerodynamic  character istics  of  the  fuselage. 

The  fuselage  aerodynamics  within  the  C81  are  represented  by 
analytic  exoressions,  which  are  third-order  functions  of  angle 
of  attack  and  sideslip.  When  these  C81  input  data  were  derived 
from  existing  wind  tunnel  map  data,  experience  showed  that  the 
C81  input  data  could  not  represent  even  small  variations  from 
smooth,  relatively  simple  curves  in  the  aerodynamic  forces  and 
moments  plotted  against  angle  of  attack  or  sideslip,  even  though 
such  variations  were  often  present  in  the  wind  tunnel  data. 
Therefore,  the  fuselage  aerodynamic  data  used  by  the  C81  were 
not  as  accurate  as  the  available  wind  tunnel  map  data. 

A C81  run  was  set  up  with  inputs  corresponding  to  conditions 
previously  used  to  model  the  CH-53A  helicopter  at  100  kt  on 
the  Sikorsky  General  Helicopter  Simulation  Program,  GENHEL. 

(The  GENHEL  program  is  described  in  References  17  and  18.) 

Table  IX  compares  the  C81  results  to  the  trim  solution  of  this 
high-speed,  nonlinear  hybrid  computer  program.  Compared  to 
the  GENHEL  solution,  the  C81  does  well  in  matching  required 
control  settings,  and  has  moderate  success  in  duplicating  the 
rotor  flapping.  The  lack  of  rotor  downwash  acting  on  the  fuse- 
lage in  the  C81  program  is  evident  again  in  the  data  of  Table  IX. 
Studying  the  values  in  this  table  of  the  rolling  and  pitching 
mo-uents  produced  by  the  main  rotor  reveals  an  inadequacy  in  the 
C81  rotor  solution.  GENHEL  correctly  predicts  the  major  moment 
contribution  about  the  pitch  axis,  while  the  C81  predicts  a 
major  contribution  about  the  roll  axis,  which  does  not  corre- 
spond (via  the  hub  moment  constant)  to  the  relative  magnitudes 
of  the  rotor  longitudinal  and  lateral  flapping  angles  prescribed 
by  the  C81.  This  result  indicates  an  inadequacy  in  the  C81  rotor 
analysis  for  predicting  the  precession  (gyroscopic)  relations  of 
the  articulated  rotor.  The  C81  rotor  analysis  was  originally 
developed  specifically  for  a teetering  rotor  system. 

The  C81  program  was  exercised  in  sideward  and  rearward  flight 
trim,  and  the  predictions  were  correlated  with  CH-53D  flight 
test  data.  The  fuselage  aerodynamic  data  were  loaded  into  the 
program  based  on  wind  tunnel  data  and  on  estimates  of  the 
helicopter  forces  and  moments  in  sideward  and  rearward  flight. 
Since  the  available  aerodynamde  data  for  the  CH-53D  in  the  wind 
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tunnel  configuration  included  the  tail,  the  tail  data  groups 
were  removed  from  the  C81  input  v/hen  running  the  sideward  and 
rearward  data  cases.  The  test  trim  data  and.  predicted  C8i 
values  for  the  CH-53D  in  sideward  flight  are  given  in  Figure  8. 
The  C81  program  showed  some  difficulty  trimming  in  sideward 
flight;  it  would  not  trim  beyond  20  kt  when  flying  to  the  left. 
The  C81  results  show  few  variations  in  roll  attitude  with  side- 
ward speed  or  direction;  nor  does  the  lateral  cyclic  control 
prediction  match  flight  test  very  accurately.  Again  the  C81 
tail  rotor  pitch  control  is  seen  to  be  low. 

Figure  9 gives  the  C81  trim  parameters  compared  to  test  data 
for  the  CH-53D  in  rearward  flight.  Evident  in  these  plots  is 
the  completely  erroneous  prediction  of  the  amount  of  lateral 
cyclic  control  required  to  trim  the  model  at  the  rearward 
speeds  studied. 

Table  X contains  the  C81  and  test  data  for  the  CH-53U  in  various 
coordinated  turns  in  the  vicinity  of  75  kt.  The  C81  coordinated 
turn  trims  were  obtained  only  after  several  attempts  to  run  the 
program  for  this  flight  condition,  due  to  difficulties  in  ob- 
taj.ning  trim  as  described  earlier.  The  C81  program  demonstrated 
reasonable  success  in  predicting  the  trim  attitudes  and  the  con- 
trol settings  in  the  turns,  considering  the  program  trim  criti- 
cisms already  described,  which  naturally  also  influence  the 
coordinated  turn  trim. 

The  overall  impression  of  the  C81  trim  analysis  as  applied  to 
articulated  rotor  helicopters  is  that  the  program  predicts  trim 
characteristics  with  moderate  success.  The  program  would 
probably  benefit  from  a more  detailed  correlation  with  articu- 
lated rotor  helicopter  trim  data  combined  with  program  refine- 
ments for  obtaining  better  correlation  in  those  areas  needing 
improvement  (such  as  the  analytic  calculations  for  tlie  aero- 
dynamic surfaces).  The  C81  appears  to  model  low -speed  charac- 
teristics well,  and  the  trim  solution  displays  proper  static 
stability  characteristics.  The  C81  trim  solution  would  be  more 
accurate  if  a more  accurate  representation  of  fuselage  aero- 
dynamics than  is  presently  available  could  be  loaded  into  the 
program.  Some  program  inadequacies  lie  in  the  rotor  solution 
for  the  articulated  system,  but  this  problem  is  not  entirely 
visible  when  studying  only  the  trim  analysis.  The  trimming 
routine,  which  brings  the  C81  to  a trim  solution,  operates  well 
in  general,  although  there  are  occasional  instances  in  which  a 
trim  solution  cannot  be  obtained  due  solely  to  some  type  of 
inefficiency  in  this  routine. 
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PERFORMANCE  CORRELATION  - RESULTS  AND  DISCUSSION 


The  C81  data  for  the  performance  correlation  were  gathered  from 
the  time-variant  portion  of  the  combined  quasi-static/time- 
variant  trim  solution  of  the  program.  Most  of  the  C81  perfor- 
mance data  shown  here  were  obtained  from  program  runs  utilizing 
four  main  rotor  mode  shapes.  Preliminary  runs  showed  little 
variation  in  power  requirements  between  cases  using  four  mode 
shapes  and  those  using  two  mode  shapes;  typically  there  was  also 
little  difference  in  power  between  the  quasi-static  and  the 
time-variant  solutions  (about  0);e  to  two  percent  difference  in 
each  instance) . Because  more  accuracy  was  expected  in  a time- 
variant  solution  using  a greater  number  of  blade  bending  modes, 
and  because  the  user's  manual  recommended  a time-variant  solu- 
tion when  doing  any  detailed  analysis  concerned  with  the  rotor 
of  an  articulated  system,  these  conditions  were  selected  when 
running  most  of  the  performance  cases.  The  use  of  five  blade 
bending  modes  was  attempted  in  the  performance  runs,  but  this 
condition  could  not  be  run  with  any  efficiency  due  to  the 
problems  discussed  in  the  previous  Trim  Correlation  section  of 
the  report.  Those  few  cases  that  did  r'on  with  five  bending 
m.odes  showed  power  requirements  of  only  one  to  two  percent 
difference  from  the  requirements  established  from  cases  using 
four  bending  modes.  The  drag  coefficient  of  the  blade  spar 
inboard  of  the  airfoil  section  was  maintained  at  the  same 
constant  value  for  all  conditions  run  in  the  performance  analy- 
sis. 

The  main  rotor  shaft  horsepower  required  in  hover  predicted  by 
the  C81  trim  analysis  for  various  CH-53A  gross  \r?eights  is 
compared  to  appropriate  test  data  in  Figure  10.  In  hover,  the 
power  requirement  predicted  by  the  quasi-static  and  by  the 
time-variant  portions  of  the  combined  quasi-static/time-variant 
C81  trim  analysis  are  identical.  The  data  of  Figure  10  show 
that  the  C81  calculates  less  power  required  to  hover  at  a given 
helicopter  gross  weight  than  is  established  by  the  flight  test 
data . 

The  main  rotor  power  coefficient  versus  tip  speed  ratio  for  the 
CH-53A  predicted  by  both  the  quasi-static  and  the  time-variant 
fjortions  of  the  C81  trim  analysis  at  forward  speed  is  compared 
to  appropriate  test  data  in  Figure  11.  Both  portions  of  the 
C81  solution  are  shown  to  provide  insight  into  the  differences 
in  performance  data  obtained  from  the  quasi-static  and  time- 
variant  solutions.  As  anticipated,  the  time-variant  analysis 
predicts  higher  power  requirements.  Both  portions  of  the  C81 
trim  solution  predict  higher  main  rotor  horsepower  than  is 
specified  by  the  test  data,  although  the  C81  data  appear  to  fit 
the  test  data  relatively  well  at  moderate  speeds. 

Figure  12  shows  the  C81  quasi-static  and  time-variant  predicted 
main  rotor  power  coefficients  (versus  main  rotor  tip  speed  ratio) 


conpared  to  test  data  for  the  S-67  helicopter.  These  data 
illustrate  that  the  C81  underpredicts  the  power  requirement  for 
the  S-67  helicopter.  As  mentioned  in  the  Trim  Correlation  sec- 
tion, the  wings  of  the  S-67  may  be  generating  too  much  lift, 
causing  the  main  rotor  power  requirement  to  be  too  low  and  thus 
incorisistent  with  the  CH-53A  comparison  in  Figure  11. 

Figure  13  repeats  the  same  C81  ouasi-static  performance  data  for 
the  CII-53A  that  is  shown  in  Figure  11  and  compares  it  to  per- 
formance data  predicted  for  the  same  conditions  by  the  Sikorsky 
Generalized  Rotor  Performance  program  (GRP) . A description  of 
the  GRP  program  can  be  found  in  References  19  and  20.  The  quasi- 
static C81  solution  was  selected  for  comparison  VN/ith  GRP  because 
the  quasi-static  analysis  more  closely  resemoles  the  GRP  analysis 
(i.e.,  a modified  single-blade  analysis)  than  does  the  time- 
variant  analysis.  Both  the  C81  and  GRP  programs  predict  the 
test  results  with  similar  accuracy;  for  the  particular  test  con- 
ditions representing  the  baseline  in  this  comparison,  the  C81 
predicts  the  moderate  speed  perform- ice  more  accurately,  while 
GRP  predicts  the  performance  at  higi.  speed  better.  The  accuracy 
of  C81  at  low  speed  is  due  to  tlie  inclusion  of  the  tip  vortex 
correction  factors  in  the  C81  calculation  of  induced  velocities 
that  are  employed  by  the  C81  analysis  even  when  tip  vortex  ef- 
fects are  not  specifically  modeled  in  the  simulation. 

Some  main  rotor  and  tail  rotor  whirl  stand  cases  were  simulated 
by  the  C81  program  by  setting  the  program  logic  switch  IPL(l) 
unequal  to  zero.  When  using  the  reduced  data  deck  for  these 
runs,  it  was  discovered  that  a tail  rotor  simulation  was  possible 
only  when  the  tail  rotor  data  were  loaded  into  the  main  rotor 
data  groups  of  the  C81.  The  C81  program  was  solved  for  various 
specified  blade  pitch  control  inputs,  obtaining  the  resulting 
rotor  thrust  values  and  the  required  rotor  power  for  the  condi- 
tions simulated.  The  main  rotor  thrust  coefficient  versus  power 
coefficient  for  the  CH-53D  predicted  by  the  C81  is  compared  to 
appropriate  whirl  test  data  in  Figure  14.  The  data  of  this 
figure  show  that  the  power  requirement  predicted  by  the  C81 
loses  accuracy  at  increasing  thrust  levels  (corresponding  to 
increasing  blade  angles  of  attack) . The  tail  rotor  thrust 
coefficient  versus  power  coefficient  for  the  CH-53A  predicted 
by  the  C81  is  compared  to  appropriate  whirl  test  data  in 
Figure  15.  The  thrust  level  for  any  given  power  predicted  by 
the  C81  for  the  tail  rotor  is  higher  than  the  measured  test 
value  of  thrust.  At  one  particular  low  setting  of  the  tail 
rotor  blade  pitch  (4  deg) , the  C81  predicted  a power  coefficient 
of  0.00038,  corresponding  to  a thrust  coefficient  of  0.0098. 

This  datum  point  was  not  consistent  vrith  the  other  thrust/power 
predictions  by  the  C81,  and  was  omitted  from  Figure  15.  The 
cause  of  this  one  inconsistent  datum  point  was  not  apparent 
upon  examination  of  the  C81  computer  run  which  produced  it. 

A number  of  C81  cases  were  run  to  iterate  the  climb  performance 
predicted  by  the  program.  At  a given  speed,  cases  were  set  up 
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requesting  trim  at  different  rates  of  climb.  Maximum  rate  of 
climb  at  a given  speed  is  based  on  available  power  and  available 
collective  pitch.  After  obtaining  trims  at  two  or  three  rates 
of  cliiab,  rate  of  climb  versus  collective  pitch  and  versus  power 
were  plotted  to  determine  the  maximum  rate  of  climb.  This  pro- 
cess v^as  then  repeated  at  a different  speed.  Figure  IG  shows 
the  maximum  rates  of  climb  obtained  in  this  manner  from  the  CBl 
at  two  altitudes  compared  to  test  data  for  tl\e  UH-5  3C.  Some 
difficulty  was  experienced  obtaining  even  a limited  number  of 
data  points  from  the  C81  because  the  program  could  not  easily 
trim  its  simulated  helicopter  in  an  attempted  climij.  The  C81 
simulated  trim  was  very  sensitive  in  the  directional  degree  of 
freedom  in  general,  including  climb  runs.  Stacked  cases,  in 
which  only  the  requested  rate  of  climb  changed,  were  also  diffi- 
cult to  run.  Often  the  quasi-static  portion  of  the  trim  could 
not  be  obtained  succescfully  for  a stacked  case,  requiring  the 
programmer  to  re-run  the  case  a number  of  times,  varying  the 
fixed  value  of  either  the  roll  or  the  yaw  attitude  until  a satis- 
factory trim  had  been  established. 

In  a (iianr.er  stinilar  ta  the  cliiilT  porforwa-.Ci^  ar^ilysis.  Aut.r^,ta- 
tive  performance  data  were  gathered  from  the  C81  for  the  HH-53C. 
At  a given  speed,  cases  were  set  up  requesting  trim  at  different 
rates  of  descent.  After  obtaining  trims  at  tv;o  or  three  rates 
of  descent,  the  rate  of  descent  varsus  the  collective  control 
was  plotted,  and  the  maximum  rate  of  descent  was  established 
based  on  the  lowest  collective  pitch  available  for  the  helicop- 
ter. (The  C81  predicted  essentially  zero  main  rotor  torque  for 
the  runs  with  low  collective  pitch.)  This  process  was  then 
repeated  at  a different  speed.  The  maximum  rate  of  descent  ver- 
sus airspeed  thus  determined  from  the  C81  is  plotted  in  Figure 
17  along  with  the  corresponding  test  results.  Unlike  the  climb 
cases,  iiO  difricuities  were  expei.lei.oeu  ii.  obtai.iiag  tiie  auto- 
rotative  trims.  The  data  of  Figure  17  show  that  the  C81  is 
highly  optimistic  in  predicting  the  autorotation  performance 
characteristic  of  the  helicopter. 

For  this  study,  the  major  prediction  effort  was  directed  within 
the  normal  flight  envelopes  of  the  subject  helicopters.  Simu- 
lated CH-53A  pull-ups  were  attempted  by  the  CBl  program  in  an 
effort  to  obtain  representative  data  points  approaching  and  on 
the  flight  boundaries,  and  the  resulting  data  points  were  to  be 
placed  on  the  empirically  determined  load  factor -velocity  (V-n) 
diagrciras  of  Reference  21.  But  such  data  were  never  obtained 
because  the  C81  program  could  not  be  trimmed  for  any  of  the 
cases  attempting  maneuver  trims.  A number  of  attempts  were  made 
to  obtain  these  maneuver  trims,  allowing  for  40  trim  iterations 
per  case.  Initial  condition  program  inputs  (particularly  con- 
trol settings)  were  varied  many  times  without  producing  a suc- 
cessful trim.  The  program  output  for  these  cases  gave  no  clear 
indication  of  v;hat  the  difficulty  was.  The  maneuver  trim  solu- 
tion difficulties  of  the  C81  program  merit  investigation  because 
this  is  a desirable  program  capability  that  would  be  frequently 
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requested  in  a heiicopter  design  effort.  If  any  further  corre- 
lation work  is  conducted  with  the  C81  program  in  the  future, 
the  maneuver  trim  data  should  be  gathered  for  comparison  to 
flight  test  data  provided  some  program  improvement  has  been  made 
allov/ing  these  cases  to  be  run. 

Tail  rotor  performance  was  not  studied  in  as  much  detail  as  main 
rotor  performance  was  studied  because  the  main  rotor  and  tail 
rotor  solutions  are  essentially  identical  in  the  C81,  and  be- 
cause tail  rotor  test  data  is  not  as  easily  accessible  as  main 
rotor  test  data.  However,  an  unusual  tail  rotor  characteristic 
was  discovered  during  the  course  of  this  correlation  effort. 
While  running  S-67  cases  on  the  C81,  the  helicopter  sometimes 
v\'ould  not  trim  when  the  tail  rotor  was  made  to  rotate  in  the 
proper  direction;  yet  when  it  was  allowed  to  rotate  in  the 
opposite  sense,  the  otherwise  identical  cases  trimmed  success- 
fully. Investigatory  cases  which  trimmed  when  either  direction 
of  tail  rotor  rotation  was  used  revealed  that  trim,  performance, 
and  response  characteristics  were  hardly  affecued  by  the  direc- 
tion of  rotation.  The  influence  of  the  direction  of  tail  rotor 
rotation  on  reaching  a trim  discouraged  any  further  detailed 
pursuit  of  tail  rotor  performance  characteristics. 

Concerning  the  general  operation  of  the  performance  analysis  of 
the  C31,  the  lack  of  rotor  wake  effects  on  the  fuselage  (as 
discussed  in  the  Trim  Correlation  section)  clearly  subtracts 
from  the  accuracy  of  the  helicopter  performance  predicted  by 
the  program.  A better  understanding  of  the  performance  analysis 
in  its  present  form  might  be  obtained  by  continued  correlation 
with  good  whirl  stand  test  data.  An  interactive  wake  routine 
would  be  needed  for  an  accurate  performance  prediction  for  a 
l.elicopter  possessing  a coaxial  or  tandem  rotor  system.  Al- 
though not  a major  criticism  of  the  C81,  the  iterative  method 
that  is  presently  required  to  determine  the  ma.ximum  rate  of 
climb  and  the  autorotative  rate  of  descent  is  cambersorne,  re- 
quiring mucii  interaction  and  plotting  of  data  by  the  program 
user.  .A  program  modification  that  would  iterate  ana  interpolate 
within  tne  program  to  determine  these  maximum  values  would  be 
useful . 


STABILITY  CORRELATION  - RESULTS  AND  DISCUSSION 


The  stability  analysis  option  provided  by  the  C81  proyram  w-^s 
used  to  gather  roots  and  stability  derivatives  for  a number  of 
Ch-53  and  S-67  correlation  cases.  The  stability  analysis  used 
by  the  C81  follows  a quasi-static  trim  of  the  case  under  inves- 
tigation. Choice  of  various  program  inputs  and  proyram  options 
(number  of  blade  bending  modes,  yawed  flow,  and  unsteady  aero- 
dynamics) were  found  to  have  essentially  no  effect  on  the  roots 
and  stability  derivatives  calculated  by  the  program.  No  diffi- 
culties related  to  the  stability  routine  of  the  program  were 
experienced  while  making  the  stability  d.  ua  runs. 

The  majority  of  the  C81  stability  runs  were  made  employing  ten 
degrees  of  freedom  for  the  helicopter:  six  body  motions,  two 

main  rotor  motions,  and  two  tail  rotor  motions.  Figures  18, 

19,  and  20  compare  tlie  roots  calculated  by  the  C81  for  the 
Cli~53A  at  three  airspeeds  (hover,  100  kt , 150  kt)  to  the  roots 
calculated  by  the  Sikorsky  GENIIEL  program  for  the  corresponding 
conditions.  The  linear  model  derived  from  the  GENHEL  analysis 
employs  nine  degrees  of  freedom:  six  body  and  three  main  rotor 
motion  (including  the  coning  degree  of  freedom) . The  roots  for 
this  model  are  shown  in  Figures  18,  19,  and  20,  and  the  six  body 
and  two  main  rotor  roots  from  the  C81  were  selected  and  shown 
for  comparison  in  these  figures.  The  roots  originating  from 
the  GENHEL  model,  in  Figures  18,  19,  and  20,  were  derived  from 
linear  derivative  models,  v/hich  were  obtained  from  GENHEL  by 
using  a system  identification  (parameter  identification)  tech- 
nique described  in  Reference  15.  Past  experience  has  shown 
that  the  body  roots  associated  with  the  nine  degrees-of -freedom 
system  identified  linear  models  are  very  similar  to  the  roots 
associated  with  linear  models  derived  by  the  familiar  perturba- 
tion technique. 

Figures  18,  19,  and  20  indicate  that  the  C81  roots  associated 
with  the  main  rotor  degrees  of  freedom  are  of  a considerably 
lower  frequency  than  the  GENHEL  rotor  roots.  The  exact  impli- 
cations of  this  characteristic  are  not  clear  in  a general 
complete-helicopter"  analysis  due  to  the  short  period  of  these 
loots,  but  the  magnitude  of  the  difference  in  the  frequencies 
of  the  rotor  roots  implies  that  some  accuracy  is  lost  in  model- 
ing the  articulated  rotor  with  the  existing  rotor  description 
used  by  the  C81. 

Figures  18,  19,  and  20  show  that  the  characteristics  of  the 
Dutch  roll  mode  and  the  phugoid  mode  predicted  by  tlie  C81  differ 
noticeably  from  those  characteristics  predicted  by  GENHEL.  The 
periods  of  the  Dutch  roll  modes  of  the  two  analyses  are  con- 
siderably different.  Figure  20  shows  that  GENHEL  models  a 
phugoid  mode  that  has  become  aperiodic  at  150  kt.  This  charac- 
teristic is  not  shown  by  the  C81  model.  The  C81  predicts  an 


aperiodic  instability  in  hover,  which  is  not  evident  in  the 
Gd'dHHL  linear  model,  v/hich  can  be  seen  .in  Figure  IB.  The  period 
and  the  time  to  double  amplitude  tor  the  phugoid  mode  ot  the 
CH-53A,  taken  from  the  data  of  Figures  13,  19,  and  20,  are  shown 
in  Figure  21  for  the  C31  and  the  GLMdEL  root  locus  solutions. 
Figures  13,  19,  and  20  illustrate  that,  in  general,  the  charac- 
teristics of  the  linear  models  derived  from  the  C81  do  not  match 
the  characteristics  of  the  linear  models  derived  from  the  GENHFL 
program. 

The  stability  derivatives  gathered  from  the  Cl'  for  the  CH-53A 
(weighing  33,500  lb  with  a 348-inch  fuselage  station  center-of- 
gravity  location)  at  100  knoti-  are  listed  in  matrix  form  in 
Table  XI,  having  been  normalized  by  mass  and  appropriate  moments 
of  inertia.  Therefore,  the  boay,  rotor  force,  ajid  moment 
variations  in  the  table  are  expressed  as  accelerations. 

Also  listed  in  Table  XI,  in  parentheses,  are  the  corresponding 
derivatives  identified  from  the  GEiJHEL  program  for  those  degrees 
of  freedom  that  are  common  to  the  two  analytic  helicopter  model 
programs.  The  derivatives  for  the  main  rotor  coning  degree  of 
freedom  derived  from  GEMHEL  are  also  shown  for  completeness. 
Appropriate  sign  convention  changes  have  been  made  to  maintain 
consistency  between  the  two  sets  of  derivatives.  The  roots 
associated  ’with  the  linear  models  described  in  this  matrix  of 
stability  derivatives  are  those  given  in  Figure  19.  The  sta- 
bility derivatives  listed  in  Table  XI  show  what  these  expressions 
look  like  for  the  ten  degrees-of-freedom  model  derived  from  the 
nonlinear  C81  helicopter  representation  and  from  the  nine  de- 
grees-of-freedora  model  derived  from  the  nonlinear  GENHEL  heli- 
copter representation.  Both  derivative  models  have  independent 
rotor  degrees  of  freedom  that  remain  constant  as  each  body  degree 
of  freedom  is  allowed  to  vary  independently.  These  linear  de- 
rivative models  are  not  comprised  of  the  stability  derivatives 
which  have  keen  more  frequently  experienced  in  the  past;  i.e., 
body  derivatives  with  respect  to  the  body  parameters,  for  which 
the  rotor  has  been  allowed  to  respond  to  and  to  retriin  before 
being  included  in  the  calculation  of  the  body  derivatives. 

Table  XII  lists  the  stability  derivatives  for  the  same  flight 
condition,  but  this  time  corresponding  to  the  six  degrees-of- 
freedom  linear  model  in  which  the  rotor  is  not  treated  as  having 
independent  degrees  of  freedom.  The  derivatives  derived  from 
the  GENIIEL  program  are  shown  in  parentheses.  These  terms  were 
determined  from  GENHEL  by  using  the  familiar  perturbation  tech- 
nique. No  means  exists  for  extracting  such  a set  of  stability 
derivatives  from  the  AGAJ74  version  of  the  C81,  so  the  C81 
derivatives  listed  in  Table  XII  were  derived  from  the  ten 
degrees-of-freedom  C81  derivatives  listed  in  Table  XI.  This 
was  done  in  the  following  manner.  Assuming  a unit  perturbation 
in  one  body  variable  (while  all  other  variable  changes  were 
fixed  at  zero) , the  four  rotor  equations  were  solved  for  four 
unknowns:  main  rotor  and  tail  rotor  fore/aft  and  lateral  flappiiag 
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angles.  Each  flapping  acceleration  and  rate  was  assumed  to  be 
zero  while  solving  the  set  of  four  equations  in  order  to  simu- 
late the  retrimmed  rotor  conditions  for  the  given  body  varial^le 
perturbation.  After  the  flapping  angles  were  solveu  for  the 
body  variable  perturbation,  these  flapping  angles  wore  multi- 
plied by  the  appropriate  body  degree-of -freedom  derivative  with 
respect  to  the  appropriate  flapping  angle.  For  each  body  degree 
of  freedom,  the  resulting  contributions  for  each  flapping  angle 
change  were  then  added  to  the  body  degree-of -freedom  derivative 
with-respect-to  the  body  variable  that  had  been  perturbed  with 
unity  magnitude.  These  sums  yield  the  stability  derivatives 
for  the  six  body  degrees  of  freedom  for  the  body  variable  being 
perturbed.  The  same  process  was  then  repeated  for  a different 
body  variable  until  the  matrix  of  Table  XII  v;as  completed. 

The  method  used  to  calculate  stability  derivatives  for  the  six 
degrees-of-f reedom  linear  model  from  tlie  C81  ten  degrees-of- 
freedom  derivatives  is  only  an  approximation  of  the  actual 
helicopter  dynamics  occurring  after  some  small  perturbation. 

The  C81  program  should  be  modified  to  provide  accurate  stability 
derivatives  for  a helicopter  with  an  articulated  rotor  v/hen  it 
is  not  desired  to  treat  the  rotor  as  having  independent  degrees 
of  freedom.  There  is  enough  confidence  in  the  concept  used  here 
to  approximate  the  C81  six  degrees-of-f reedom  stability  deriva- 
tives to  permit  these  data  to  be  compared  to  the  corresponding 
derivatives  derived  from  the  GEhilEL  program,  in  an  effort  to 
determine  what  portions  of  the  C81  description  of  an  articulated 
helicopter  need  further  investigation  or  improvement. 

The  important  stability  derivatives  for  describing  helicopter 
stability  characteristics  modeled  by  the  C81  were  selected  from 
Table  XII  and  were  compared  to  corresponding  derivatives  obtained 
from  C,ENHEL.  Derivatives  that  showed  more  than  small  differences 
from  each  other  were  studied  to  determine  what  might  be  respon- 
sible for  these  differences.  Insight  gained  from  all  otlier  C81 
correlation  efforts  (trim,  performance,  time  histories,  and  rotor 
loads)  were  applied  in  this  stability  study,  and  references  to 
the  higher  degrees-of-freedom  derivatives  of  Table  XI  were  made 
whenever  these  more  specific  (but  less  familiar)  derivatives 
could  assist  in  the  analysis. 

The  lack  of  a routine  for  calculating  the  influence  of  rotor 
downwash  on  the  fuselage,  and  the  accuracy  lost  by  representing 
the  fuselage  aerodynamics  by  an  analytic  expression  rather  than 
by  wind  tunnel  map  data  have  some  influence  on  the  discrepancies 
seen  in  the  C81  calculation  of  X-force  with  vertical  velocity 
and  Z-force  with  longitudinal  velocity  derivatives  of  Tables  XI 
and  XII.  The  pitching  moment  with  vertical  velocity  is  also  a 
function  of  the  aerodynamic  characteristics  of  the  fuselage. 

The  orientation  of  the  rotor  thrust  relative  to  the  fuselage 
affects  all  of  the  major  longitudinal  stability  derivatives 
which  the  C81  does  not  predict  accurately.  All  those  stability 
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derivatives  that  are  a function  of  rotor  thrust  orientation 
appear  to  be  predicted  with  little  accuracy  by  the  C81.  The 
difference  in  signs  betv/een  the  C81  and  GENHEL  derivatives  for 
rotor  pitching  moment  with  lateral  flapping  and  rot  rolling 
moment  with  longitudinal  flapping  (see  Table  XI)  shows  an 
inaccuracy  in  rotor  cross-coupling  in  the  C81.  The  C81  trim 
solution  of  the  rotor  and  the  cross-coupling  effects  due  to 
rotor  precession  calculated  by  the  C81  program  for  an  articu- 
lated rotor  system  lead  to  enough  questions  to  merit  a detailed 
investigation  of  all  contributions  to  these  calculations.  Such 
an  investigation  should  include  a comparison  of  the  various 
compon:;nts  of  each  rotor  force  and  moment  to  the  corresponding 
components  calculated  by  some  established  baseline  analysis  for 
articulatfid  rotors  (such  as  the  GENHEL  program)  . This  type  of 
detailed  ..nves tigation  would  probably  require  entering  the  C81 
program  and  breaking  the  solution  at  various  points  in  the 
programmed  equations.  The  technique  used  to  approximate  the 
articulated  rotor  solution  within  the  C81  should  be  further 
investigated  if  the  C81  rotor  solution  in  the  program  is  sub- 
jected to  detail  ^ examination. 

The  damping  and  control  power  expressions  from  the  C81  in 
Table  XII  compare  well  with  those  calculated  by  GENHEL  except 
in  the  directional  degree  of  freedom.  The  yawing  moment  with 
yaw  rate  derivative  has  the  wrong  sign  (opposite  direction) , 
while  the  yawing  moment  with  tail  rotor  pitch  is  too  high  as 
predicted  by  the  C81.  The  side  force  and  the  rolling  moment 
derivatives  with  respect  to  tail  rotor  pitch  are  also  too  high. 
The  tail  rotor  modeled  by  the  C81  program  therefore  must  produce 
too  much  change  in  thrust  in  response  to  a change  in  tail  rotor 
pitch.  Other  major  lateral-directional  derivatives  that  are  not 
predicted  accurately  by  the  C81  are  the  rolling  moment  with  yaw 
rate,  and  the  yawing  moment  with  lateral  velocity  and  with  roll 
rate.  These  stability  derivatives  are  functions  of  the  change 
in  tail  rotor  thrust  with  sideslip  and  of  the  change  in  side 
force  produced  by  the  vertical  tail  with  sideslip. 

Some  C81  stability  data  runs  were  made  in  order  to  collect 
CH-53D  and  S-67  root  locus  data  for  comparison  to  roots  esti- 
mated from  flight  test  data.  Figures  22  and  23  show  the  body 
roots  of  the  CH-53A  at  100  kt  and  150  kt  obtained  by  system- 
identifying  flight  test  data.  The  roots  obtained  from  the  C81 
program  runs  simulating  these  same  flight  conditions  also  appear 
in  the  figures.  Roots  gathered  by  semiempirical  metliods  show 
some  difference  from  analytically  derived  roots.  If  the  GENHEL 
set  of  roots  from  Figures  19  and  20  is  compared  to  the  two  sets 
of  roots  in  Figures  22  and  23,  respectively  (the  gross  weight 
and  center-of -gravity  locations  of  the  flight  conditions  repre- 
sented by  these  sets  of  data  are  slightly  different) , it  is  seen 
that  the  GENHEL  roots  more  closely  resemble  the  test  data  roots 
than  do  the  C81  roots. 


Few  root  locus  data  are  available  for  the  S-67  helicopter. 

Roots  derived  from  test  data  are  shown  in  Figures  24,  25,  and 
26  in  order  to  provide  some  such  data  for  the  S-67  in  this 
reports  These  roots  were  obtained  by  studying  time  history 
test  history  test  data  and  measuring  the  frequency  and  period 
of  reponses  to  pulse  input  disturbances.  The  frequency  and 
period  information  was  transferred  to  root  locus  plots.  Roots 
gathered  from  C81  data  runs  corresponding  to  the  same  flight 
conditions  are  also  included  in  Figures  24,  25,  and  26. 

The  stability  characteristics  of  the  C81  linear  representations 
of  the  CH-53  and  S-67  helicopters,  as  interpreted  from  the  roots 
and  derivatives  calculated  by  the  C81  program,  can  only  reflect 
the  characteristics  of  the  nonlinear  model  from  which  the  sta- 
bility derivatives  and  roots  were  derived.  The  roots  and  deriva- 
tives cannot  be  expected  to  be  highly  accurate  if  some  portions 
of  the  nonlinear  helicopter  representation  need  improvement. 

The  roots  and  derivatives  are  used  to  indicate  the  areas  of 
weakness  (and  strength)  in  the  model  from  which  they  were  de- 
rived. The  mechanics  of  the  stability  analysis  within  the  C81 
for  gathering  roots  and  derivatives  works  well  and  relatively 
quickly.  Besides  a needed  option  for  obtaining  the  stability 
derivatives  of  an  articulated  rotor  helicopter  selected  without 
independent  rotor  degrees  of  freedom,  the  ability  to  obtain 
stability  derivatives  at  any  time  during  a time-variant  maneuver 
solution  might  also  prove  to  be  a useful  addition  to  the  C81 
program. 
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TIME  HISTORY  CORRELATION 


- RESULTS  AND  DISCUSSION 


Maneuver  cases  for  the  H-53  and  S-67  helicopters  were  run  to 
demonstrate  the  response  characteristics  of  the  C81  modeled 
helicopter  via  the  resulting  time  history  data.  Pulse,  step, 
and  other  miscellaneous  control  inputs  from  reliable  test  data, 
representing  a variety  of  selected  aircraft  configurations  and 
flight  conditions,  were  reproduced  as  accurately  as  possible  and 
submitted  as  inputs  into  the  C81  program  during  Its  maneuver 
mode  of  operation.  All  helicopter  maneuver  responses  simulated 
by  the  C81  and  compared  to  flight  test  data  were  for  basic  heli-- 
copter  responses  without  stability  augumenta tion.  Experience 
showed  that  the  number  of  main  rotor  bending  modes  used  in  the 
maneuver  cases  had  negligible  effect  on  response  time  hi  stories, - 
all  of  the  time  history  data  runs  referenced  in  the  figures  were 
made  using  two  main  rotor  normal  modes  (hinged  rigid  flapping 
and  lagging) . Most  of  the  data  cases  were  run  using  the  time- 
variant  solution  of  the  C81,  although  the  responses  differed 
very  little  from  the  responses  obtained  from  the  quasi-static 
maneuver  solution.  A few  time  history  data  cases  were  run  using 
the  quasi-static  solution  in  order  to  study  the  time  history 
characteristics  of  the  rotor  flapping  angles.  Yawed  flow  was 
employed  in  all  the  time  history  data  runs,  although  unsteady 
aerodynamics  were  requested  in  only  a few  of  the  runs.  The 
effects  of  yawed  flow  and  unsteady  aerodynamics  on  the  param- 
eters under  study  were  negligible. 

Quasi-static  maneuvers  for  the  CH-53A  at  approximately  33,000 
lb  and  aft  center-of-gravity  location  were  run  on  the  C81  pro- 
gram at  hover,  100  kt,  and  150  kt  to  compare  the  resulting  heli- 
copter responses  to  flight  test  data.  In  particular,  the  rotor 
longitudinal  and  lateral  flapping  angle  responses  following 
longitudinal  stick  pulse  inputs  were  to  be  compared;  these  angles 
are  available  only  from  the  quasi-static  solution  of  the  C81. 

Few  of  these  flapping  response  correlation  data  are  shown  in  the 
report,  The  problem  with  these  (and  other)  maneuver  time  history 
C81  runs  is  that  immediately  upon  beginning  the  time  history 
solution,  nonzero  angular  accelerations  are  calculated  by  the 
program.  Often  these  accelerations  are  of  large  magnitudes  (up 
to  70  deg/sec ^ for  one  of  the  HH-53C  cases) . These  erroneous 
inputs  appear  as  discontinuous  disturbances  to  the  helicopter. 
Some  cases  show  the  helicopter  responding  in  the  wrong  direction; 
reaction  to  the  erroneous  acceleration  input  completely  overrides 
the  control  stick  input.  These  initial  nonzero  accelerations 
occur  frequently  for  both  quasi-static  and  time-variant  time 
history  solutions.  These  acceleration  disturbances  usually 
occur  at  the  very  beginning  of  the  time  history  solution,  al- 
though in  isolated  instances  they  appeared  later  in  the  time 
history.  For  all  the  time  history  data  cases  that  were  run  in 
this  study,  the  helicopter  was  quasi-static  trimmed  well  within 
a very  tight  tolerance,  usually  within  10  pounds  or  foot-pounds. 


For  this  reason,  and  because  the  problem  occurred  for  quasi- 
static maneuvers  as  well  as  tiiae-varian t maneuvers,  the  nonzero 
accelerations  are  most  likely  caused  by  a program  memory  location 
error . 

The  quasi-static  longitudinal  pulse  cases,  which  were  run  on  the 
C81  for  the  33  ,000  lb  CH-53A,  typically  had  very  poor  pitch  atti- 
tude and  load  factor  responses  in  comparison  to  the  test  results. 
Figure  27  presents  data  from  the  C81  run  with  the  least  inaccu- 
rate aircraft  responses,  from  among  these  CH-53A  cases.  Due  to 
the  inexact  duplication  of  the  helicopter  response  by  the  C81 
program  for  this  case,  firm  conclusions  cannot  be  drawn  on  the 
rotor  flapping  response  predicted  by  the  C81.  However,  the 
resulting  longitudinal  flapping  r -ponse  predicted  by  the  C81 
for  this  case  peaked  at  2.5  deg  v spared  to  about  1 deg  from 
test  data),  while  lateral  flapping  was  cross-coupled  by  about 
1 deg  in  the  prediction  (compared  to  almost  zero  response  in  the 
test  data) , as  seen  in  Figure  27. 

Noticeably  absent  from  the  C81  program  output  are  values  of  the 
flapping  and  lagging  angles  of  individual  rotor  blades  during 
the  time  history.  Even  if  these  angles  were  only  approximated 
by  the  blade  and  rotor  head  geometry,  based  on  the  generalized 
coordinates  of  the  closest  blade  segment  outboard  of  the  flapping 
and  lagging  hinge  locations,  these  data  would  prove  useful  when 
using  the  C81  to  model  helicopters  with  articulated  rotor  systems 
and  should  be  made  a^'-ailable  from  the  C81. 

Figures  28  through  36  show  the  time  history  responses  predicted 
by  the  C81,  compared  to  test  data,  for  pulse  and  step  control 
inputs  to  the  HH-53C  helicopter.  These  data  are  shown  here 
even  though  initial  nonzero  angular  accelerations  occasionally 
occur  in  the  C81  runs.  All  C81  accelerations  at  time  zero  in 
tliese  figures  are  shown  as  zero  because  the  program  prints  these 
zero  values  when  beginning  a maneuver.  But  the  large  nonzero 
value  for  acceleration  seen  at  the  next  data  point  printed  in 
the  computer  output  (requested  as  early  as  0.004  sec  into  the 
maneuver  for  some  of  the  data  runs  made)  indicates  that  an 
erroneous  input  is  immediately  entering  the  maneuver  solution. 

Figure  28  , 29  , and  30  illustrate  the  actual  and  the  C81  predicted 
responses  of  the  HH-53C  at  31,000  lb  in  hover  to  longitudinal, 
lateral,  and  directional  pulse  control  inputs,  respectively. 
Figure  30  illustrate  the  effects  oi  the  large  erroneous  accelera- 
tions sustained  in  roll  and  yaw  at  the  beginning  of  this  C81  run. 
The  quasi-static/ time-variant  trim  condition  solutions  generated 
by  the  C81  prior  to  initiating  the  maneuvers  illustrated  by  the 
time  histories  of  Figures  29  and  30  were  identical;  yet  the  non- 
zero acceleration  disturbance  that  occurred  during  the  lateral 
input  case  was  negligible  in  comparison  to  that  for  the  pedal 
input  case.  Irregularities  in  the  time  histories  of  Figure  29 
near  4 seconds  elapsed  time  indicate  that  some  extraneous 


acceleration  distu^Lanco  can  occur  at  some  time  other  than  zero 
in  the  C81  maneuver  solution. 
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The  unusual  behavior  of  the  pitch  acceleration,  both  near  time 
zero  and  near  three  seconds,  corresponding  to  a forv/ard  stick 
pulse  input  predicted  by  the  C81,  as  seen  in  Figure  20,  prevents 
any  conclusions  frorti  being  drawn  from  the  data  for  this  case. 
However,  a divergence  in  yaw  can  be  seen,  and  the  directional 
degree  of  freedom  may  be  highly  coupled  to  the  roll  response. 

The  data  of  Figure  20  indicate  decent  correlation  between  C31 
and  test  data  for  a right  stick  pulse  input  except  that  the 
roll  attitude  response  corresponding  to  the  changes  in  roll  rate 
and  acceleration  predicted  within  the  program  was  unusual.  The 
data  of  Figure  30  does  not  lead  to  any  useful  information  about 
the  C81  maneuver  solution  because  of  the  severe  initial  inputs 
in  roll  and  yaw  accelerations. 

At  31,000  lb  in  hover,  the  actual  and  the  C81  predicted  responses 
of  the  HH-53C  to  longitudinal,  lateral,  and  directional  step 
control  inputs  are  shovm  in  Fiaures  31,  32,  and  33,  respectively. 
Despite  an  initial  roll  acceleration , which  was  far  too  large, 
the  response  of  the  helicopter  predicted  l^y  the  C81  for  an  aft 
stick  step  input,  shown  in  Figure  31,  closely  resembled  the 
actual  test  response.  The  variation  in  pitch  responses  may  have 
been  due  to  a minor  but  erroneous  biasing  of  pitch  acceleration 
at  time  zero.  In  this  example  the  coupled  roll  response  appears 
to  be  proper,  once  the  bias  due  to  different  trim  values  of  roll 
attitude  is  considered.  The  data  of  Fioure  32  show  that  the 
roll  rate  and  acceleration  predicted  by  the  C81  were  too  large 
within  the  2 seconds  following  a ricjht  stick  step  input.  Figure 
30  also  indicates  that  the  coupled  pitch  response  to  the  roll 
input  was  not  in  the  proper  direction.  Not  much  can  be  concluded 
from  the  data  of  Figure  33  about  the  response  to  a pedal  step 
predicted  by  the  COl  because  of  the  large  initial  disturbance 
in  yaw  acceleration. 

Figures  34,  35,  and  36  show  the  actual  and  the  CBl  predicted 
response  of  the  HH-53C  at  41,000  Ih  at  113  kt  to  longitudinal, 
lateral,  and  directional  step  control  inputs,  respectively. 

The  pitch  acceleration  predicted  by  the  CBl  for  a forward  stick 
step  input,  as  shown  in  Figure  34,  appears  to  have  been  delayed 
compared  to  the  flight  test  results.  The  change  in  load  factor 
predicted  by  the  CBl  v-zas  not  as  great  as  the  change  in  the 
actual  data.  Figure  35,  illustrating  the  response  to  a left 
stick  step,  also  indicates  a delay  in  the  roll  response  pre- 
dicted by  the  C81  when  compared  to  the  test  data.  The  C81  data 
of  Figure.  36,  for  a right  pedal  step  input,  was  contaminated  by 
roll  and  yaw  acceleration  disturbances  at  about  0.5  second 
elapsed  tine,  which  occurred  for  no  apparent  reason. 

The  Col  data  runs,  which  yield  the  data  of  Figures  28  throuch 
36,  employed  a 15-doaree  change  in  main  rotor  azimuth  between 


successive  program  calculations  during  each  time  history  solu- 
tion. The  C81  data  plotted  in  these  figures  were  drawn  from 
the  computer  printouts,  which  requested  listings  of  the  maneuver 
solution  only  after  every  60-degree  change  in  rotor  azimuth. 

Some  CALCOMP  plots  were  requested  using  a more  frequent  data 
sampling  interval  (every  30  degrees  in  azimuth) . The  CALCOMP 
plots  occasionally  revealed  a 6/rev  oscillatory  content  in 
angular  accelerations  during  the  maneuvers.  This  oscillatory 
content  was  not  exhibited  in  the  computer  printed  data  because 
of  the  infrequent  sampling  of  data. 

Some  S-67  time  histories  solved  by  the  C81  were  gathered  and 
compared  to  test  data,  but  these  runs  also  showed  high  nonzero 
accelerations  at  the  beginning  of  the  computed  maneuver  solution. 
A symmetrical  dive  and  pull-up  at  180  kt  was  run  on  the  C81. 
Comparing  the  program  data  to  test  data  over  25  seconds  of  time 
history  response  showed  the  C81  to  predict  excessive  pitch 
angles  (as  much  as  60  deg  nose  down) . The  normal  load  factor 
and  the  change  in  airspeed  solved  by  the  C81  were  also  extremely 
inaccurate.  About  30  seconds  of  response  following  an  aft  cyclic 
stick  pulse  input  solved  by  the  C81  were  studied  for  the  S-67  at 
140  kt.  These  data  showed  roll  to  be  highly  coupled  in  the  solu- 
tion, reaching  excessive  roll  angles  (60  deg)  in  comparison  to 
test  data.  Here,  too,  the  load  factor  and  airspeed  responses 
were  erroneous.  The  pitch  response  was  divergent,  showing  no 
oscillatory  content  as  seen  in  the  test  data.  The  root  locus 
solution  from  the  C81  for  exactly  the  same  aircraft  configuration 
and  flight  condition  (Figure  25)  predicts  long-term  longitudinal 
stability  characteristics  different  from  those  shown  by  the  non- 
linear C81  time  history  solution  of  this  data  case.  Another 
S-67  maneuver,  which  was  simulated  by  the  C81  at  140  kt , was  a 
pedal  doublet.  Studying  about  20  seconds  of  response  showed  the 
yaw  response  of  the  helicopter  solved  by  the  C81  to  be  excessive 
(40  deg)  and  divergent  (unlike  test  data) . The  roll  response  of 
the  helicopter  was  highly  coupled  and  excessive  (35  deg) . The 
time  history  of  the  sideslip  predicted  by  the  C81  did  compare 
well  with  test  results.  It  should  be  remembered  that  the  initial 
accelerations  of  the  S-67  for  all  the  cases  discussed  here  were 
excessive,  and  tiiis  may  have  played  a major  part  in  the  ensuing 
responses. 

Although  most  of  the  C81  time  histories  gathered  in  this  study 
were  apparently  contaminated  by  the  unusual  acceleration  solu- 
tions, some  data  were  presented  and  discussed  in  this  report  any- 
how, when  it  was  felt  that  the  data  might  indicate  characteris- 
tics of  the  basic  C81  maneuver  solution.  Some  of  the  time  his- 
tory data  suggests  that  the  C31  program  may  be  moderately  suc- 
cessful in  predicting  articulated  rotor  helicopter  maneuver  re- 
sponse characteristics.  (However,  the  computer  time  required  for 
time  history  solutions  is  excessive  and  is  discussed  in  a later 
section  of  this  report.)  Pitch-roll  coupling  and  yaw  sensitivity 
app>ear  to  be  areas  needing  improvement.  Detailed  quantitative 
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analysis  of  the  C81  maneuver  solution  should  be  postponed  until 
the  excessive  accelerations  experienced  in  the  solution  are 
corrected  or  explained- 
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ROTOR  LOADS  CORRELATION  - RESULTS  AND  DISCUSSION 


Data  for  the  C81  rotor  loads  correlation  wore  collected  by 
running  C81  maneuver  cases,  with  no  disturhino  inputs,  for  a 
duration  equivalent  to  one  rotor  revolution,  and  printouts  of 
blade  moment  solutions  were  requested  after  trim  and  periodically 
durincf  the  rotor  revolution.  A considerable  effort  v/as  expended 
running  introductory  and  comnarative  cases  to  determine  how  to 
best  set  up  the  cases  to  prodaico  tlie  most  accurate  rotor  loads 
solutions  from  the  COl  for  an  articulated  rotor  helicopter. 

After  establishinn  criteria  for  selecting  input  data  and  pronram 
options,  data  cases  were  run,  and  C81  rotor  loads  data  v;ere 
collected  for  comparison  to  test  and  other  analytic  data.  7>,11 
rotor  loads  cases  run  in  this  study  were  for  the  CH-53A  heli- 
copter . 

Comparing  data  from  a number  of  CHI  runs  indicated  that  the 
! blade  moment  distribution  solution  of  this  pronram  is  highly 

dependent  upon  (and  far  too  sensitive  to)  the  degree  of  acr  acy 
used  in  expressing  the  blade  normal  mode  shapes  and  mass  di.  , ri- 
bution.  Figure  37  shows  the  distribution  of  the  blade  bear.wise 
! 1 bending  moment  (at  60  deg  rotor  azimuth  location)  calculated  by 

■ the  C81  for  a numljor  of  cases.  Table  -XIII  presents  the  blade 

I normal  modes  shapes  and  mass  distributions  used  for  these  cases. 

! Tv/o  mode  shapes  (rigid  flapping  and  lagging)  were  used  in  each 

case  for  the  CH-53A  at  32,550  lb  in  hover.  (The  rigid  lagging 
. mode  shapes  wore  identical  to  the  rigid  flapping  moc3e  shapes 

listed  in  Table  XIII.)  Case  A in  the  figure  illustrates  the 
1 solution  for  an  ideal  blade  with  exactly  linear  mode  shapes 

(calculate  ' to  four  significant  figures)  and  uniform  mass  distri- 
I bution.  Case  B is  the  solution  for  a blade  which  also  has  uni- 

( form  mass  distribution,  but  whose  mode  shape  generalized  coordi- 

i nates  are  only  approximately  linear.  (The  inexactness  of  the 

mode  shape  at  5%  radius  for  case  B is  the  result  of  an  attempt 
^ ) to  accurately  represent  the  location  of  the  flapping  hinge. 

jj  This  is  another  problem  which  will  be  discussed  shortly.)  Com- 

^ paring  the  data  of  cases  A and  B in  Figure  37  shows  that  spikes 

i appear  in  the  moment  solution  wherever  the  mode  shape  is  not 

^ exact,  and  that  the  moment  distribution  across  the  blade  can  be 

'1  biased  by  these  discontinuities.  Case  C illustrates  the  moment 

■ 'I  distribution  when  exactly  linear  bending  mode  shapes  are  used, 

‘ but  when  the  mass  distribution  for  the  blade  is  nonuniform. 

Case  C is  compared  to  case  A for  this  purpose. 

IVhen  the  C81  program  is  used  to  model  a helicopter  with  an 
articulated  rotor,  the  locations  of  flapping  and  lagging  hinges 
are  specified  through  the  normal  mode  shape  input  data  for  the 
rotor  blades.  The  flapping  hinge  offset  from  the  center  of 
* rotation  of  the  rotor  for  the  CH-53A  is  5.56%  of  the  rotor 

radius.  Due  to  the  restriction  in  C81  input  and  program  equa- 
tions regarding  the  fixed  number  of  blade  segments  (20) , each 
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of  equal  length  (5%  radius)  , the  exact  offset  for  an  articulated 
rotor  may  be  loaded  into  the  program  only  if  the  offset  is  an 
integer  multiple  of  5%  radius.  Cases  b and  D of  Figure  37 
illustrate  unsuccessful  attempts  to  accurately  locate  the  true 
hinge  offset  by  loading  mode  shapes  into  the  program  which  re- 
flect the  exact  location  of  this  hinge.  These  cases  were 
unsuccessful  because  the  C81  cannot  interpolate  the  proper  off- 
set location  between  blade  station  mode  shape  input  data.  Case 
D represents  the  moment  solution  when  the  blade  input  data 
represent  the  exact  nonuniform  mass  distribution  and  the  pre- 
cise offset  hinge  location  of  the  CH-53A  main  rotor  blade. 

Another  problem  pertaining  to  the  flapping  hinge,  discovered 
while  making  these  and  other  investigatory  C81  rotor  loads  runs, 
involves  the  C81  solution  of  the  bending  moment  at  tlie  flapping 
hinge.  This  moment  should  always  be  zero.  The  only  condition 
for  which  the  C81  would  solve  a zero  moment  at  the  hinge  (within 
some  computational  error  tolerance)  was  for  the  hinge  to  be 
located  exactly  at  a blade  segment  radial  station,  with  exactly 
linear  blade  mode  shapes,  and  with  uniform  blade  mass  distribu- 
tion (case  A in  Figure  37).  All  other  cases  that  were  run  (many 
more  than  are  shown  in  Figure  37)  could  produce  no  better  an 
average  value  of  hinge  moment  over  one  revolution  than  1000  to 
2000  'n-lb. 

Investigatory  C81  rotor  loads  solution  runs,  which  were  made  in 
addition  to  those  included  in  Figure  37  showed  that  the  number 
of  normal  mode  shapes  used  in  the  runs  and  the  forward  speed 
of  the  helicopter  do  not  affect  the  results  seen  in  Figure  37 
or  discussed  in  the  previous  paragraphs.  Properties  similar  to 
those  discussed  were  also  exhibited  for  fhe  chordwise  moment 
distribution,  as  shown  in  Figure  38. 

The  majority  of  the  rotor  load  runs  made  in  this  study  used 
the  shaft-free  solution  of  the  C81,  wherein  the  entire  helicopter 
is  trimmed.  But  the  shaft-fixed  solution  was  also  investigated, 
though  only  in  a few  computer  runs.  Here  the  shaft  angle  is 
specified  and  only  the  rotor  is  trimmed,  either  for  specified 
flapping  angles,  or  for  specified  control  angles.  Table  XIV 
presents  the  rotor  trims  obtained  froia  various  shaft-fixed  solu- 
tions, compared  to  free  aircraft  (shaft-free)  solutions,  and  to 
test  data.  These  C81  solutions  correlated  well  with  test  data 
on  longitudinal  flapping  and  control  angles,  although  the  lateral 
flapping  and  control  angles  predicted  by  the  C81  were  not  always 
as  accurate.  The  shaft-fixed  solutions  were  obtained  using  two 
normal  modes  (flapping  and  lagging).  A quasi-static  trim  could 
not  be  reached  when  shaft-fixed  trims  were  attempted  using  five 
blade  bending  modes. 

Runs  were  made  using  two  and  five  normal  mode  shapes  to  compare 
the  resulting  blade  bending  moment  distribution  along  the  blade 
radius.  The  steady  beamwise  moment  at  the  flapping  hinge  was 
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considerably  smaller  when  two  modes  were  used  instead  of  five 
(-1G00  in-lb  compared  to  8800  in-lb  for  the  CH-53A  at  34,120  lb 
and  159  kt) ; but  the  oscillatory  content  of  bending  moments  was 
much  too  small  everywhere  else  along  the  blade  when  two  modes 
were  used.  The  vibratory  bending  moments  solved  by  the  C81  when 
five  nodes  were  used  matched  test  data  much  more  closely  than 
they  did  when  two  modes  were  used.  Therefore,  it  was  decided 
to  use  five  mode  shapes  in  the  rotor  loads  data  runs  being 
ccnf'ared  to  test  and  other  analytic  data. 

The  C81  data  runs  made  for  the  purpose  of  collec'ina  blade 
bending  moments  for  comparison  to  test  and  analytic  results  v/ere 
shaft-free  solutions  using  five  main  rotor  blade  bending  modes. 
Control  system  deformation,  not  explicitly  available  in  the  C31 
program,  was  represented  in  the  torsional  bending  mode  shape 
loaded  into  the  runs.  Yawed  flow  and  unsteady  aerodynamics  v/ere 
included  in  tliese  runs.  The  flapping  hinge  was  located  at  5% 
radius  rather  than  5.5G%  in  order  to  avoid  biasing  the  moment 
dis  ti.  i)jution , which  would  have  occurred  had  an  attempt  been  made 
to  accurately  locate  the  flapping  hinge.  The  lagging  hinge  was 
located  at  a designated  blade  station  for  the  same  reason.  Mode 
shapes  were  loaded  into  these  runs  as  accurately  as  was  possible, 
and  the  correct  (nonuniform)  blade  mass  distribution  was  used. 

Three  C81  rotor-load  data  runs  were  compared  to  test  results 
for  the  CH-53A  helicopter.  These  were  for  34,120  l.b  in  hover, 
33,070  lb  at  159  kt,  and  40,580  lb  at  137  kt  (all  at  a neutral 
center  of  gravity).  These  conditions  cover  low  and  high  weight, 
as  well  as  hover  and  high  speed.  The  C81  data  gathered  at  137 
kt  and  159  kt  were  also  compared  to  analytic  solutions  from  the 
Sikorsky  Y200  normal  modes  blade  aeroelastic  analysis  computer 
program.  (The  Y200  pronram  is  described  in  References  22  and 
23.)  The  data  from  the  Y200  runs  v/ere  previously  obtained  i.i 
conjunction  with  the  work  of  Reference  24.  (Due  to  the  solution 
technique  used  in  program  Y200,  the  rotor  thrust  obtained  in  it=: 
solutions  arc  about  6%  higher  than  the  thrust  values  for  the 
test  and  the  C81  conditions.) 

Figure  39  compares  the  C81  and  Y200  solutions  of  the  steady 
contribution  of  beamwise  bending  moment  versus  blade  radius  to 
test  data  for  the  CH-53A  at  159  kt.  The  C81  more  closely  pre- 
dicts the  maximum  steady  beamwise  bonding  moment.  The  Y200 
predicts  the  tip-down  steady  moment,  whereas  the  C81  does  not. 
Figure  40  is  similar  to  Figure  39  except  that  it  compares  the 
vibratory  content  of  the  beamwise  bending  moment.  Here,  the 
Y200  predicts  the  proper  trend  in  the  radial  distribution  of 
vibratory  beamwise  bending  moment,  while  the  C81  does  not. 

Figure  41  compares  the  harmonic  content  of  the  vibratory  beam- 
wise  bending  moment  near  the  radial  station  where  this  moment 
is  highest.  (The  test  and  Y200  data  actually  correspond  to  79^ 
radius  but  are  compared  to  the  CBl  values  at  30?,  since  the  C81 
program  fixes,  by  5?  increments,  the  radial  stations  at  which 


bondinq  noncnt  data  will  bo  available.)  The  two  proorans  pre- 
dict lov;- f roquency  (1/rov  and  2/rev)  beamwiso  bondinq  noment 
anplitutd.es  with  conparable  accuraev.  The  Y20n  program  predicts 
the  approximate  amplitudes  of  the  high-frequency  beamwise  bend- 
ing moment  components.  The  C81  underpredicts  these  vibratory 
components.  Figure  42  compares  time  )iistories  of  tlie  the  total 
beamwise  moment  over  ono  rotor  revolution  at  the  critical  blade 
station  (00%  radius)  from  the  CBl,  the  Y2nn,  and  the  flight  test 
data.  This  tir.e  historv  is  predicted  accurately  by  the  Y200 
analysis.  The  corresponding  COl  time  history  does  not  exhibit 
the  proper  phase  relationshiyi  or  frequency  content. 

The  steady  chordwise  bending  moments  (versus  blade  radius)  pre- 
dicted by  the  COl  and  the  Y200  analyses  are  compared  to  the  test 
d.ata  in  Fig'ire  43  for  the  CII-53A  at  159  kt.  The  CBl  approximates 
the  radial  distribution  of  the  steady  chordwise  bending  moment 
slightly  better  than  doo.s  the  Y200.  Figure  44  compares  vibratory 
chordwise  bending  moments  solved  by  the  CBl  and  the  Y200  to  test 
data.  Both  the  CBl  and  the  Y200  analyses  underprediict  tlie  ampl  i- 
tude  of  chordwise  vibratory  bending  moment.  The  CBl  yields  a 
slightly  better  prediction  of  flight  test  results.  Taking  a 
closer  look  at  the  vibratory  component  at  the  critical  blade 
station,  the  harmonic  content  of  the  vibratory  chordwise  bending 
moment  at  55%  radius  predicted  by  the  CBl  and  the  Y200  is  com- 
pared to  test  results  in  Figure  45.  The  test  data  in  this  figure 
show  the  major  fourth  harmonic  contribution,  which  is  a charac- 
teristic of  this  CH-53A  rotor  blade.  Although  the  Y200  under- 
predicts the  magnitudes  of  the  harmonic  content,  the  Y200  analy- 
sis does  indicate  that  the  fourth  harmonic  is  the  major  contri- 
butor to  this  vibratory  bending  moment.  The  CBl  results  do  not 
predict  this  characteristic,  but  it  should  be  noted  here  that 
Y200  employs  a detailed  variable  inflow  model  (including  inter- 
active calculation  of  vortex  influence  on  the  inflow) , while  the 
CBl  uses  a much  simpler  analytic  expression  of  inflow.  The  time 
.histories  of  the  total  chordwise  bending  moment  predicted  by  the 
(.  81  and  the  Y200  during  one  rotor  revolution  at  55%  radius  are 
com.pared  to  test  data  in  Figure  4fi.  As  with  beamwise  m.oment, 
the  Y200  program  more  accurately  predicts  the  time  history  of 
blade  chordwise  bending  moment  than  does  the  CBl  program. 

A limited  number  of  blade  torsiona]  moment  data  were  available 
for  the  cases  studied.  Figures  47  and  48  show  the  steady  and 
vibratory  content,  respectively,  of  the  torsional  moments  pre- 
dicted by  the  CBl  and  the  Y200,  compared  to  CH-53A  test  data, 
for  the  159  kt  case.  Figure  49  compares  test  results  to  the  CBl 
and  the  Y200  solutions  for  the  harmonic  content  of  vibratory 
torsional  moment  near  15%  radius.  Both  programs  had  low  predic- 
tions of  the  high-f rec(uency  (near  G/rev)  components  of  torsional 
moment  associated  with  retreating  blade  stall  flutter. 

Similar  to  the  rotor  loads  comparison  just  discussed  for  the 
CH-53A  at  33,070  lb  and  159  kt,  CBl  and  Y200  analyses  were  also 
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run  and  compared  to  test  data  for  the  CH-53A  at  137  kt,  .t  at 
a higher  rotor  loading  corresponding  to  40,630  lb  gross  weight. 
With  few  exceptions,  the  C81  and  the  Y200  each  performed  in  the 
same  manner  when  predictina  blade  bending  moments  for  the  higher 
weight  case  as  they  did  for  the  lower  weinht  predictions  illus- 
trated in  Figures  3'^  through  4^.  The  differences  in  the  results 
for  the  two  conditions  were  that  for  the  hiaher  weiaht  case,  the 
C81  and  the  Y200  predicted  almost  identical  steady  chordv;ise 
bending  moment  distribution  with  span  (both  underpredicting  the 
flight  test  results) . The  vibratory  chordwise  bending  moments 
versus  span  predicted  by  these  two  analyses  also  were  almost 
identical  (and  once  again  both  underpredicted  corresponding  test 
results).  For  the  high  weight  case,  the  vibratory  torsional 
moment  distribution  with  span  predicted  by  the  Y200  was  slightly 
higher  than  the  prediction  of  the  C81.  Both  the  C81  and  the 
Y200  pr-  grams  underpredicted  the  magnitude  of  vibratory  torsional 
moments  for  the  high  thrust  condition. 

The  C81  rotor  loads  solution  for  the  C1!-53A  in  hover  (at  34,120 
lb)  was  compared  to  test  data  in  the  same  manner  as  was  done  at 
159  kt.  Tliis  comparison  led  to  basically  the  same  results  and 
conclusions  as  were  drawn  from  the  comparison  made  at  159  kt. 

Few  differences  were  observed  between  these  two  cases.  In  hover, 
the  C81  underpredicted  the  steady  beamv;ise  bending  moment  near 
the  blade  tip  and  predicted  far  too  low  vibratory  beamwise 
bending  moments  over  the  entire  blade  span.  For  this  case 
the  C81  mistakenly  saw  essentially  only  a first  harmonic  contri- 
bution to  the  vibratory  beamwise  bending  moment,  and  grossly 
underpredicted  the  magnitude.  In  hover,  the  C81  underpredicted 
the  steady  chordwise  bending  moments  over  most  of  the  blade, 
except  at  the  tip  where  it  overpredicted.  For  this  hover  case 
the  C81  underpredicted  the  vibratory  chordwise  bending  moments 
over  the  entire  span  of  the  blade  and  underpredicted  all  the 
harmonic  components  of  this  moment  at  the  critical  blade  section. 

To  summarize  the  evaluation  of  the  C81  program  for  predicting 
blade  bending  moments  for  an  articulated  rotor  helicopter,  the 
program  solution  is  overly  sensitive  to  variations  in  normal 
mode  s)iapes  and  the  mass  distribution  of  the  blade.  The  articu- 
lated hinges  cannot  be  accurately  located  due  to  the  fixed  seg- 
ment length  restriction.  The  C81  program  solved  an  erroneous 
nonzero  moment  solution  at  the  hinge,  and  this  characteristic  of 
the  solution  creates  some  lack  of  confidence  in  bending  moment 
solution  elsewhere  along  the  blade.  In  hover,  the  C81  underpre- 
dicts the  vibratory  content  of  blade  bending  moments.  At  high 
speed,  with  higher  bending  modes  inputted  to  the  program  in  addi- 
tion to  the  two  rigid  flapping  and  lagging  modes,  the  C81  gener- 
ally predicts  the  proper  magnitude  of  the  vibratory  content  of 
bending  moments,  while  predicting  the  steady  contribution  with 
only  moderate  success.  The  C81  does  not  predict  the  radial  dis 
tribution  of  moment  well  where  this  distribution  shows  local 
peaks  or  sudden  changes  in  value.  Using  the  analytic  expression 


within  the  C31  program  to  model  the  rotor  inflow,  the  program 
does  not  accurately  recognize  the  proper  harmonic  content  of  the 
vibratory  bending  moments.  Time  histories  of  the  blade  bending 
moments  do  not  show  the  proper  frequency  content  and  are  not 
correctly  synchronized  with  the  rotor  azimuth  location  of  the 
blade . 

.■Regarding  the  rotor  loads  input/output  options  of  the  C31,  more 
freedom  in  locating  blade  stations  at  which  the  rotor  loads  can 
be  analyzed  is  aesirable.  In  addition  to  blade  bending  moments, 
other  variables  that  are  involved  in  the  rotor  loads  analysis 
(such  as  blade  stresses)  should  be  laade  available  to  the  user  if 
these  variables  are  accessible  or  computable  within  tlie  existing 
program.  The  user  is  presently  restricted  in  the  number  of  times 
he  may  request  the  printout  option  which  reveals  the  bending 
moments  along  each  radial  station  of  each  blade  during  a time 
history  solution;  this  option  should  be  improved  to  allow  the 
user  to  see  these  values  as  frequently  as  he  desires  during  a 
rotor  revolution  following  the  trim  solution. 
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DISCUSSJLON  J3F  J2™ERAj;^_PRO_rnmM_  CHARACTERISTICS 

Thi55  section  of  the  report  discusses  those  C81  input,  opera- 
tional, output,  and  auxiliary  program  characteristics,  that 
proved  to  be  cumbersome  or  were  found  to  need  improvement.  The 
topics  discussed  here  were  not  included  in  previous  sections  of 
the  report  either  because  their  detailed  discussion  was  not 
necessary  in  explaining  the  immediate  topics  or  because  these 
topics  apply  to  the  CSl  program  in  general. 

The  time  required  bv  the  central  processing  unit  (CPU)  of  the 
computer  to  execute  the  major  analyses  of  the  CBl  program  is  an 
important  characteristic  of  the  program.  The  required  CPU  time 
for  a solution  determines  the  cost  of  the  solution  and  therefore 
is  a very  practical  consideration  when  evaluating  the  usefulness 
of  the  C81  program  as  an  articulated  rotor  helicopter  design 
tool.  High  CPU  requirements  by  the  CSl  are  rticularly  unac- 
ceptable if  other  analyses  are  available  thu  >rovide  solutions 
that  can  be  obtained  using  less  CPU  time  than  CBl.  The  version 
of  tlie  CSl  programmed  at  Sikorslcy  using  the  HX  compiler  ran 
twice  as  fast  as  the  ver.sion  compiled  with  the  C2  compiler. 

Virtual  storage  of  the  CSl  program  within  the  computer  was  used 
in  preference  over  real  storage  so  that  the  computer  could  be 
shared  with  other  programs  while  running  the  CSl.  Typical  quasi- 
static trims  and  stability  analyses  renuired  about  one  minute  of 
CPU  time.  This  v.’as  found  to  be  an  acceptable  renuirement. 
Ouasi-static/time-variant  combined  trims,  hov/ever,  required  from 
four  to  eight  minutes  of  CPU  tine.  The  required  CPU  time  deoends 
greatly  on  the  rotor  azimuth  increment  used  in  consecutive  solu- 
tions v/ithin  a time-variant  trim  or  during  any  maneuver  time  his- 
tory solution.  CPU  renuirements  were  found  to  depend  only 
slightly  on  other  program,  trim  routine  inputs  varied  during  the 
course  of  this  study.  Maneuver  tim.e  histories  were  found  to 
require  excessive  amounts  of  CPU  time.  Four  seconds  of  time 
history  required  about  25  CPU  minutes  when  the  rotor  azimuth 
increment  of  15  degrees  v/as  used,  while  about  50  CPU  minutes 
v/ere  required  for  the  same  time  history  when  the  rotor  azimuth 
increment  used  was  9 degrees.  In  comparison,  the  Sikorskv  CENHEL 
program  runs  fast  enough  to  be  used  in  a real-time  sim.ulation 
of  helicopter  maneuver  responses.  The  time  needed  for  the  CSl 
time  history  solution  did  not  depend  greatly  on  whether  a quasi- 
static maneuver  or  a tim.e-variant  maneuver  was  requested.  A typi- 
cal rotor  loads  analysis  run  made  in  this  study,  including  a 
maneuver  time  history  solution  over  one  rotor  revolution  (with 
a 5-deg  rotor  azimuth  increment) , required  about  20  minutes  of 
CPU  tim.e.  This  requirement  was  somevfhat  high  in  com.parison  to 
the  CPU  requirements  of  the  Sikorsky  Y200  program  for  similar 
solutions . 

Certain  articulated  rotor  helicopter  characteristics  are  not 
represented  in  the  CSl  program.  The  CSl  model  has  no  rotor 


coniny  ^.jyroo  of  froedoi.i  anc  tiiorefore  calculates  no  coniny 
anjle.  There  is  no  ^irovision  in  the  J81  to  induce  a blade 
lay- to-pi tcJi  coupliny  ainy.:.  deyardiny  more  yeneral  helicOi>tsr 
and  rotorcraft  inputs,  tije  C81  has  no  input  tnat  expresses  tlie 
uynaniic  pressure  loss  at  the  tail  rotor.  The  stability  auycienta- 
tion  system  and  the  control  iuixiny  offered  in  the  coiitrol  sys- 
tem input  yroup  encompasses  many  possible  conf  iyura  tions ; iiow- 
ever,  they  are  not  flexi'de  enouyh  to  allow  .lany  comraon , 
straiyhtforward  stability  and  mixing  inputs  to  uj  represented  in 
the  program.  An  improvement  is  needed  her  that  would  allow  the 
easy  input  of  any  stability  aug.nen tation  or  control  laixing  ar- 
rangeraent  to  the  C81  rotorcraft  analytic  model. 

dome  of  tlie  necessary  input  data  to  tiie  Cdl  proved  to  be  difii- 
cult  to  prepare  or  were  cumbersome  or  confusing  to  use.  The 
fuselage  aerodynamic  data  were  converted  from  wind  tunnel  map 
uata  into  input  cojasistent  with  the  Col  program  by  using  the 
\S812.\  program.  This  program  was  written  in  PL/1  language  and 
was  supplied  to  the  contractor  by  the  Lustis  directorate-  The 
AS312A  was  developed  to  convert  fuselage  aerodynamic  data  into 
input  consistent  witdi  the  C81.  The  7\S812A  program  was  difficult 
to  run  successfully.  The  restrictive  input  data  requirements  of 
the  AS312A  laade  it  necessary  to  load  in  more  aerodynamic  input 
data  than  were  often  available.  The  engineer  v/as  forcec  to 
estimate  some  bivariate  aerodynamic  map  data  to  satisfy  the  in 
put  requirements  of  i\S0127w  This  auxiliary  program  should  have 
aiore  flexible  input  requirements. 

Another  point  pertaining  to  tlie  fuselage  aerodynamic  data  input 
of  the  Cdl  concerns  the  format  of  these  uata  witn  regard  to  tlie 
nominal  and  high  angle  of  attack  equations.  It  is  unclear  what 
input  is  required  for  different  flight  conditions,  particularly 
v/hen  also  considering  rearward  or  sideward  flight  conditions. 

A change  in  xnput  format  or  an  improved  description  in  the 
user's  manual  of  the  present  format  would  be  hel^gful. 

ilegarding  the  normal  modes  input,  the  blade  bending  mode  shape 
data  inputtee  to  the  C81  first  had  to  be  derived  from  another 
source.  Sikorsky  program  Y178  (coupled  noriual  modes)  >vas  used 
to  solve  the  normal  mode  shapes.  It  would  be  convenient  if  the 
laode  shapes  were  calculated  within  the  Col,  so  that  the  appro- 
priate input  data  could  be  loaded  into  the  C31  program  rather 
than  into  some  auxiliary  program. 

Son\e  errors  and  inconveniences  are  evident  in  the  C81  program 
printout  and  manuals.  The  sign  of  the  tail  rotor  thrust  printed 
xn  the  trim  output  never  changes,  even  when  tiie  tail  rotor  is 
rotated  180  deg  from  one  side  of  the  tail  boom  to  the  other. 

Tne  printout  frequently  does  not  include  titles,  descriptions, 
or  units  and  thus  requires  tiie  program  user  to  consult  the 
manuals  much  too  often.  The  units  of  some  of  the  printed  data 
(particularly  the  stability  derivatives)  change  from  place  to 
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place  in  the  printout.  This  situation  v;ould  be  confusing  even 
if  the  units  were  listed  in  the  printout. 

A particularly  annoying  and  confusing  inconvenience  in  the  pro- 
gram regards  the  different  conventions  used  in  numbering  the 
blade  stations  appearina  in  the  manuals,  the  printout,  and  the 
CALCOMP  plots.  Sometimes  the  blade  station  reference  begins  at 
the  root  of  the  blade,  sometimes  at  the  tip.  References  to 
these  stations  appear  in  many  places  in  the  program  and  should 
he  made  consistent.  One  example  is  that  the  printer  plot  lists 
the  wrong  blade  station  for  the  plotted  rotor  bending  moment 
data  calculated.  Another  examp"!  ' is  the  procram  input  XMR(15), 
which  actually  specifies  the  bl  uie  station  at  which  bending 
moments  are  calculated  when  se  inent  numbering  begins  at  the  tip, 
rather  than  at  the  root  as  e 'lained  in  the  user's  manual. 

Another  inconvenience  of  the  C81  printout  is  the  omission  of 
rotor  moment  summaries  in  the  force  and  moment  summary,  printed 
out  while  the  trim  iterations  are  being  performed.  Such  infor- 
mation can  prove  to  be  helpful  in  analyzing  a case  for  which 
the  user  is  having  difficulty  obtaining  a trim. 

The  engineer's  manual  could  be  more  informative  if  it  included 
descriptions  of  the  entire  C81  program,  rather  than  just  dis- 
cussing those  changes  made  since  the  program  was  last  documented. 
This  would  eliminate  the  need  for  the  curious  user  t--  consult 
a number  of  referenced  reports,  which  are  often  not  readily 
available . 

The  plotting  routine  included  in  the  C81  program,  which  is 
interactive  with  the  CALCOMP  system,  operated  successfully  and 
was  interfaced  with  the  CALCOMP  hardware  and  software  systems 
with  little  difficulty.  The  CALCOMP  plots,  which  were  obtained 
while  running  some  of  the  rime  history  data  cases  for  this 
study,  proved  to  be  useful  working  plots  and  aided  in  the 
correlation  effort.  The  CALCOMP  plots  as  generated  by  the  C81 
routine,  however,  are  not  of  report  quality,  and  could  be  im- 
proved. For  instance,  when  the  user  requested  more  than  one 
variable  to  be  plotted  on  the  same  graph,  it  v/as  often  difficult 
to  determine  which  plot  corresponded  to  v/hich  variable.  Also, 
when  it  was  necessary  to  use  more  than  one  scale  for  the  graph 
when  plotting  more  than  one  variable,  the  proper  scale  pertain- 
ing to  each  variable  was  not  obvious.  After  printing  the  scale 
for  the  plotted  variable,  the  pen  was  sometimes  not  commanded 
to  be  lifted  from  the  paper  before  beginning  the  time  history 
trace,  thus  producing  an  unwanted  line  on  the  graph.  The  pro- 
gram user  does  not  have  control  of  the  time  scale  selected  for 
the  CALCOMP  plots.  Giving  the  user  this  option  could  frequently 
prove  to  be  useful  while  gathering  desired  CALCOMP  plots. 

In  concluding  this  discussion  of  the  general  program  operating 
characteristics,  anyone  interested  in  using  the  C81  program 


should  be  aware  of  the  fact  that  the  program  is  large  and 
detailed,  and  considerable  tiir.e  is  required  to  prepare  input 
and  to  become  familiar  witli  the  operation  of  the  program.  The 
associated  manuals  for  the  C31  program,  standing  alone,  do  pro- 
vide enough  information  to  allow  a new  user  to  set  up  and  run 
the  program  v'ithout  any  other  assistance  or  reference,  altliough 
close  coordination  with  the  developer  or  major  users  of  the  pro 
gram  (such  the  Eustis  Directorate) 
rec^uired  by  the  new  user  in  obtaining 


greatly  reduces  the  effort 
successful  program  runs. 


CONCLUSIONS 


It  was  found  that  it  was  possible  to  nake  the  C81  proyrani 
operate  and  converge  to  trim  solutions  for  an  articulated  rotor 
helicopter  and,  within  certain  limitations,  provide  answers 
comparable  to  those  of  other  analytical  methods.  i'ne  major 
limitations  identified  in  this  study  of  the  aGAJ74  version  of 
the  C81  were: 

1.  Various  C31  analyses  (particularly  tLae-variant  triins 
and  time  history  solutions)  require  exceptionally  iiiyli 
computer  run  time.  I'he  run  times  cannot  bo  substantially 
reduced,  even  when  requesting  minimal  amounts  of  tne 
available  sophistication  and  complexity  of  the  program 
analyses . 

2.  The  C81  program  tri-rmsing  routine  is  not  efficient  enough 
to  produce  successful  trims  as  frequently  (i.e.,  few 
iterations)  as  desired  in  a user-oriented  helicopter 
design  tool  program. 

3.  The  general  trim  characteristics  of  articulated  rotor 
helicopters  in  forward  flight  (with  the  exception  of  roll 
characteristics)  are  predicted  with  moderate  success  by 
the  C81.  Tlie  program’,  does  not  accurately  predict  siuov.’ard 
and  rearward  flight  trii,.s. 

4.  The  performance  characteristics  predicted  by  the  C31  do 
not  correlate  as  well  as  desireu  with  test  data  for  some 
conditions  specified  within  the  operational  flight  envelope. 
The  required  i>iain  rotor  collective  pitch  is  not  predicted 
accurately  in  these  circorastances . The  progra.a  has  diffi- 
culty solving  maneuver  trir.-.s  (pull-ups,  etc.)  near  flight 
envelope  boundaries. 

5.  for  the  articulated  rotor  lielicopter  cases  studied  in  this 
correlation  effort,  the  stability  charac ter istics  predicted 
by  the  Cdl  prograr.;  were  inaccurate  and  generally  unaccep- 
table . 

G.  The  aynaniic  response  ciiaracteris tics  luodeled  by  the  CBl 
iielicopter  siaiulation  show  excessive  pitcn-roll  coupling 
and  yaw  sensitivity. 

7.  discontinuous  pulse  or  step  inputs  in  aelicopter  accelera- 
tion frerjuently  occur  during  the  solution  of  tlie  d81 
maneuver  analysis,  probably  originating  in  tne  i,,atnciv.atics 
of  the  prograidTiing  technique  rather  than  from  tno  governing 
dyna;aic  analysis,  thereby  contaminating  tii-e  history 
response  data. 


d.  Jladc  I’cndiny  lUoruents  solved  by  the  Cbl  rotor  loads 

analysis  are  too  sensitive  to  i:>lade  i..ass  distribution  ani.i 
.'•’.ode  shai^jo  input  data. 

9.  In  tile  present  fon.i  of  the  program,  tiie  articulated  hinpes 
of  tile  rotor  cannot  be  accurately  located  in  the  031 
iielicopter  .lodel.  Oiie  existing  techni'^ue  for  solving  rotor 
i.ioi.ients  witiiin  the  C81  does  not  solve  tiie  proper  bending 
moment  at  the  articulated  iiinges. 

10.  'i'iie  magnitudes  of  tiie  steaay  and  vil^ratory  bending  moments 
outboard  of  the  articulated  hinges  solved  by  tiie  031  rotor 
loads  analysis  are  comnaralile  to  soD.utions  o.jtuined  from 
otiier  rotor  loads  analyses.  uov/ever , tiie  time  iiistory 
solutions  of  these  bending  uioi.ients  about  tjie  rotor  azir.utli 
inaccurately  predict  when  peak  values  occur  anu  inaccura ue ly 
describe  tlie  freguoncy  content  of  tiie  luonients. 
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RECOMMENDATIONS 


The  incorporation  of  the  C81  program  as  an  articulated 
rotor  helicopter  design  tool  over  existing  separate 
helicopter  performance,  handling  qualities,  and  rotor 
loads  analyses  is  not  recommended.  This  recommendation 
is  made  in  light  of  the  comparison  of  articulated  rotor 
helicopter  characteristics  as  predicted  by  the  C81  to  test 
data  and  data  accumulated  from  other  analyses. 

If  the  C81  program  is  to  be  maintained  and  improved  as  a 
general  rotorcraft  global  program,  it  is  recommended  that 
this  program  undergo  a further  detailed  comparison  with 
test  data,  provided  program  improvements  have  been 
implemented  that  correct  the  problems  uncovered  in  this 
study.  The  general  trim  characteristics  should  be  checked 
again,  and  characteristics  established  in  maneuver  trims 
near  helicopter  operating  boundaries  should  be  correlated 
with  flight  test  data.  Performance  with  respect  to  high- 
speed power  requirements  and  autorotative  capabilities 
should  be  reevaluated.  More  accuracy  in  stability  deriva- 
tives and  root  locus  prediction  is  required.  Time  history 
response  sliould  be  studied  again,  checking  for  helicopter 
cross-coupling  effects  and  overall  sensitivity  to  control 
inputs.  Rotor  moments  and  cross-coupling  should  be 
rechecked,  as  well  as  the  time  history  solution  of  blade 
bending  moments.  To  operate  as  an  efficient  and  truly 
modular  global  program,  it  is  also  recommended  that  the 
running  tim.es  of  each  of  the  various  CRl  analyses  be 
reduced,  particularly  when  the  program  options  that  are 
requested  do  not  reouire  high  degrees  of  available  program 
sophistication . 
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TABLE  VIII.  - Continued 


TABLE  IX. 


Trim  Parameters 

deg 
0^,  deg 
deg 

MR  Thrust,  lb 
MR  X-force,  lb 
MR  side  force,  lb 
MR  Roll  moment,  ft-lb 
MR  Pitch  moment,  ft-lb 
MR  Torque,  ft-lb 


-0.11 


-0.07 


■K-nH-t 


TABLE  XI.  ROTOR  AND  BODY  STABILITY  C 

OBTAINED  FROM  THE  C81  AND 


Vari  abl  e 

u 

V 

w 

P 

q 

r 

0 

a 

Degree 

freedom 

ft/sec 

ft/sec 

ft/sec 

rad/sec 

rad/sec 

rad/sec 

rad 

rad 

u,  ft/sec"^ 

-.049 

.017 

.14 

1 .27 

-10.8 

-1.12 

0 

-32.1 

_ 

(-.005) 

( . 0004 ) 

( .094) 

(-.126) 

(-9.97) 

( .865) 

(0) 

(-32.2) 

( 

0 

V,  ft/sec' 

.003 

-.182 

-.06 

9.85 

.194 

-174. 

32.1 

0 

(.016) 

(-.057) 

(.0022) 

(6.63) 

(-.659) 

(-167.) 

(32.2) 

( .054) 

( 

0 

w,  ft/sec^' 

.151 

-.05 

-1.21 

-7.55 

168. 

4.3 

.607 

-2.3 

(-.124) 

(-.007) 

(-.744) 

(-6.12) 

(193.) 

(4.4) 

(.89) 

(-1.44) 

1 

rad/sec^ 

.0006 

-.034 

-.015 

-.56 

-.15 

1.30 

0 

0 

(-.001) 

(-.016) 

(-.010) 

(-•172) 

(.017) 

(.443) 

(0) 

(0) 

0 

q,  rad/sec" 

.0014 

-.0014 

.0003 

-.054 

-.294 

.049 

r» 

u 

0 

(-.0005) 

(.0023) 

(.0016) 

(-.006) 

(-.121) 

(.0045) 

(0) 

(0) 

0 

r,  rad/sec^ 

-.0012 

.025 

-.006 

.177 

.172 

-1  .43 

0 

0 

(-.0016) 

(.008) 

(-.0002) 

(■111) 

(.06) 

(-.728) 

(0) 

(0) 

«*1mR’ 

-.042 

.181 
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( -.038) 

(8.46) 
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.0925 

.471 

27.4 
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-.717 
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0 

0 
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0 

0 
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O0 

^OMR>  rad/sec2 
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(-.432) 
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(0) 
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Stability  derivatives  in  parentheses  are  from  GENHEL  analysis. 
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ROTOR  AND  BODY  STABILITY  DERIVATIVES  FOR  THE  CH-S3A  AT 
OBTAINED  FROM  THE  C81  AND  THE  GENHEL  ANALYSES 
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Second  derivatives  (in  parentneses)  are  from  GENHEL  analyst 


CH-53A  MAIN  ROTOR  BLAUE  TLAPPING  MODE 
SHAPES  AND  MASS  DISTRIBUTION  USEU  IN 
INVESTIGATORY  C8l  ROTOR  LOADS  SOLUTION  RUNS 
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Figure  1.  CH-53A  Articulated-Rotor  Helicopter. 
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Figure  7.  Concluded 
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Figure  9. 


C81  Rearward  Flight  Trim  va  Airspeed,  Compared 
With  Flight  Teat  Data  for  the  CH-53D  at  46,000 
Lb,  Forward  FSCG  » 328.3  (104%  h^  ° 2000  Ft) 
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Figure  12.  C81  Level  Flight  Power  Requirement  va  A( 

Ratio,  Compared  With  Flight  Test  Data  f( 
S-67  at  16,500  Lb,  102.5%  Nr , h<j  “ 0 Ft 
(FSCG  o 272.3,  Cw  “ 0.00505). 
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Figure  14,  C81  Prediction  of  Main  Rotor  Thrust  Coeffi' ient 

vs  Power  Coefficient,  Compared  with  Whirl  Stand 
Test  Data  for  the  CH-53D  Main  Rotor  Out  of  Ground 
Effect,  at  100%  Nr,  and  at  Standard  Sea  Level 


2400 


AIRSPEED  , V , KNOTS 


C81  Maximum  Rate  of  Climb  Performance  vs  Airspe^ 
Compared  With  Flight  Test  Data  for  the  HH-53C  a 
35,000  Lb  at  Various  Density  Altitudes  (100% 
FSCG  » 340) . 


3200 


2400 


---  C8I 

O FLIGHT  TEST 


AIRSPEED  , V , KNOTS 


Figure  17 


C81  Maximum  Rate  of  Descent  in  Autorotation  va 
Airspeed,  Compared  With  Flight  Test  Data  for  the 
HH-53C  at  35,000  Lb,  100%  Ny,  h^  » 4500  Ft 
(FSCG  o 340). 
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Figuro  18.  C81  Main  Rotor  and  Body  Root  Locus,  Compared  With 

GEINHEL  Analysis  Data  for  the  CH-53A  in  Hover 
(33,500  Lb,  FSCG  = 348). 
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Figure  22. 


C81  Body  Root  Locus  Compared  With  flight  T^^t 
System-Identified  Data  for  the  CH-53A  at  100  Kt 
(35,000  Lb,  FSCG  “ 352). 
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Figure  24.  C81  Body  Root  Locus  Compared  With  Flight  Test 

Time-Hiatory-Fitted  Data  for  the  S-67  at 
100  Kt  (14,800  Lb,  FSCG  = 276). 
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Figure  25, 


C8l  Body  Root  LOCOB  Co^»«^  S-67®at  Uo\t 

TiiBB-History-Pltted  oata  tor  the  s. 

(14,800  Lb,  FSCG  « 276). 
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Figure  26.  C81  Body  Root  Locus  Compared  With  Flight  Teat 

Time-History-Fitted  Data  for  the  S-67  at  180  Kt 
(14,800  Lb,  FSCG  = 276) . 
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p-  Figure  27.  C81  Maneuver  Response  Time  History  Comparison 

« With  Flight  Test  Data  for  the  CH-53A  at  100  Kt, 

Following  a Forward  Stick  Pulse  Input  (33,539 
Lb,  FSCG  = 347.7) . 
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C81  Maneuver  Response  Time  History  Comparison 
With  Flight  Test  Data  for  the  HH-53C  in  Hover, 
Following  a Forward  Stick  Pulse  Input  (31,000 
Lb,  FSCG  = 328) . 
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Figure  30.  C81  Maneuver  Response  Time  History  Comparison 

With  Flight  Teat  Data  for  the  HH-53C  in  Hover, 
Following  a Left  Pedal  Pulse  Input  ^"^1,000  Lb, 
FSCG  “ 328)  . 
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Figure  31.  C81  Maneuver  Response  Time  History  Comparison 

With  Flight  Test  Data  for  the  HH--53C  in  Hover  ^ 
Following  an  Aft  Stick  Step  Input  (31,000  Lb, 
FSCG  =»  328)  . 


123 


C8l 


FLIGHT  TEST 


AO  -T 


I 

-40 


ELAPSED  TIME  , SECONDS 


C81  Maneuver  Response  Time  History  Comparison 
With  Flight  Test  Data  for  the  HH-53C  in  Hover, 
Following  a Right  Lateral  Stick  Step  Input 
(31 ,000  Lb,  FSCG  « 328) . 
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Figure  32.  Continued. 
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C81  Maneuver  Response  Time  History  Comparison 
With  Flight  Teat  Data  for  the  KH-53C  In  Hover 
Following  a Right  Pedal  Step  Input  (31,000  Lb 
FSCG  « 328) . 
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Figure  34.  C81  Maneuver  Response  Time  History  Comparison 

With  Flight  Test  Data  for  the  HH-53C  at  113  Kt, 
Following  a Forward  Stick  Step  Input  (41,000  Lb 
FSCG  » 328) . 
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Figure  34,  Concluded. 
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Figure  35.  C81  Maneuver  Response  Time  History  Comparison 

With  Flight  Test  Data  for  the  HH-53C  at  113  Kt, 
Following  a Left  Lateral  Stick  Step  Input 
(41,000  Lb,  FSCG  = 328)  . 
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Figure  36.  C81  Maneuver  Response  Time  History  Comparison 

With  Flight  Teat  Data  for  the  HH-53C  at  113  Kt 
Following  a Right  Pedal  Step  Input  (41,000  Lb, 
FSCG  = 328) . 
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Y200  analysis 


Blade  Beeunwise  Bending  Moment  at  Approximately 
80%  Blade  Span,  Compared  With  Flight  Test  and 
Y200  Analysis  Data  for  the  CH-53A  at  159  Kt, 

33,070  Lb  (103%  KV/  = 3000  Ft,  FSCG  » 336,  Level  Flight) 
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f Figure  43.  C81  Rotor  Loads  Solution  of  Blade  Steady  Chordwise 

Bending  Moment  vs  Span,  Compared  With  Flight  Test 
and  Y200  Analysis  Data  for  tiie  CH-53A  at  159  Kt, 
33,070  Lb  (103%  Nj.,  h^  = 3000  Ft,  FSCG  «=  336, 

Level  Flight). 
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Figure  47.  C81  Rotor  Loads  Solution  of  Blade  Steady 

Torsional  Moment  vs  Span,  Compared  With  Flight 
Test  and  Y200  Analysis  Data  for  the  CH-53A  at 
159  Kt,  33,070  Lb  (103%  N , h^  “ 3000  Ft, 

FSCG  “ 336,  Level  Flight). 
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* Figure  48.  C81  Rotcr  Loads  Solution  of  Blade  Vibratory 

Torsional  Moment  vs  Span,  Compared  With  Flight 
Test  and  Y200  Analysis  Data  for  the  CH-53A  at 
159  Kt,  33,070  Lb  (103%  Nj.,  h^j  = 3000  Ft, 

FSCG  » 336,  Level  Flight). 
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rotor  coning  angle,  deg 

main  rotor  lateral  cyclic  control,  deg 

rotor  longitudinal  flapping  angle,  deg 

main  rotor  longitudinal  cyclic  control,  deg 

rotor  lateral  flapping  angle,  deg 

buttTine,  in. 

center-of -gravity 

rotor  power  coefficient 

rotor  thrust  coefficient 

rotor  loading  (weight)  coefficient 

stationline  (fuselage  station),  in, 

fuselage  station  center-of-gravity 
location  relative  to  a fixed  reference,  in. 

density  altitude,  ft 

main  rotor 

actual  rotor  speed  referred  to  design 
rotor  speed,  pet 

normal  load  factor  at  helicopter  center 
of  gravity,  g's 

notation  for  normal  load  factor  in 
velocity  - load  factor  plots,  g's 

helicopter  roll  angular  velocity,  deg/sec 

torque,  ft-lb 

helicopter  pitch  angular  velocity,  deg/sec 
helicopter  yaw  angular  velocity,  deg/sec 
tail  rotor 

helicopter  longitudinal  velocity,  ft/sec 
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helicopter  total  velocity,  kt 
helicopter  lateral  velocity,  ft/sec 
waterline,  in. 

helicopter  normal  velocity,  ft/sec 

helicopter  longitudinal  body  axis 

helicopter  lateral  body  axis 

helicopter  vertical  body  axis 

fuselage  angle  of  attack,  deg 

helicopter  body  pitch  attitude,  deg 

main  rotor  collective  blade  pitch  control,  deg 

tail  rotor  blade  pitch  control,  deg 

advance  ratio:  helicopter  velocity  over 

main  rotor  tip  speed 

helicopter  body  roll  attitude,  deg 

helicopter  body  yaw  attitude,  deg 

main  rotor  azimuth  location  of 
reference  blade,  deg 

denotes  first  derivative  with  respect 
to  time 

denotes  second  derivative  with  respect 
to  time 

denotes  blade  pitch  angle  at  75%  rotor  radius 

denotes  a parameter  referenced  to  the  helicopter 
body  axis  system 

denotes  a main  rotor  parameter 

denotes  a tail  rotor  parameter 

denotes  a change  from  the  trim  value  of  a 
parameter 
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